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Because of their small electromagnetic corrections, the isospin-breaking decays ry — >■ Stt seem to be 
good candidates for extracting isospin-breaking parameters ~ (m^ — mu). This task is unfortunately 
complicated by large chiral corrections and the discrepancy between the experimentally measured 
values of the Dalitz parameters describing the energy dependence of the amplitudes of these decays 
and those predicted from chiral perturbation theory. We present two methods based on an analytic 
dispersive representation that use the information from the NNLO chiral result and the one from 
the measurement of the charged ?7 — )■ Stt decay by KLOE together in a harmonized way in order 
to determine the value of the quark mass ratio R. Our final result is R = 37.7 ± 2.2. This value 
supplemented by values of nis/fh or even m and rus from other methods (as sum-rules or lattice) 
enables us to obtain further quark mass characteristics. For instance the recent lattice value for 
majrh ~ 27.5 leads to Q = 23.1 ± 0.7. We also quote the corresponding values of the current masses 
m„ and nid- 



I. INTRODUCTION 



The masses of the light quarks are fundamental free parameters of the standard model. Since quarks are confined 
inside hadrons, there is no direct method for their measurement. The only method of determination is a comparison of 
the theoretical prediction for some observable that depends on the quark masses with the corresponding experimental 
value. For that end we need a framework, in which the quark masses occur explicitly, and which can make predictions 
for such observables with sufficient precision. Because of quark confinement and the fact that these masses are very 
small in comparison to the typical hadron scales, perturbative quantum chromo dynamics (QCD) cannot play such a 
role and we need to employ a non-perturbative method. Nowadays, the natural candidates for such approaches are 
lattice QCD [1] and chiral perturbation theory (ChPT) 

While the mass and the isospin average mass of and , which is denoted as 

m= , (1) 

have become accessible through the recent lattice simulations (among others 041]) and agree well with the independent 
determination via QCD sum rules ^9|-(TT|, the extraction of the individual masses m„ and rrid from the lattice is still 
polluted by various simplifications of the electromagnetic effects that have to be made in isospin-breaking simulations 
(cf. e.g. |12|). Therefore, if we want to determine the individual masses niu and m^, at the moment ChPT seems to 
be the more promising approach. 

The most suitable processes for studies of isospin breaking within ChPT in the mesonic sector are 77 — > Stt decays. 
Since these decays violate G-parityQ they have to proceed via isospin-breaking effects. There are two mechanisms 
of this breaking, either through the electromagnetic (EM) interactions, which are proportional to the electric charge 
squared, or through the isospin-breaking mass difference between u and d quarks, 

nifcBi^) = {d{x)d{x) - u{x)u{x)) . (2) 

Even though the EM interactions have a sizable effect on the difference m^± — to^o and on the pion decay constant 
Ftt, it has turned out that their influence on the t] ^ Stt decay amplitudes is very small |13l - [l5j . Hence, "Hq^q 
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^ Or equivalently, the decay would bo forbidden as a result of isospin conservation and charge conjugation invariance (C-invariance). Indeed, 
the final state has to have the total isospin / = and is therefore totally antisymmetric with respect to the permutation of the three 
pions (the only allowed state is then Tr+vr^Tr'^). Due to Bose symmetry, the corresponding amplitude is then totally antisymmetric under 
exchanges of the momenta of these three pions. On the other hand, according to C-invariance, the amplitude has to be symmetric with 
respect to the exchange of the momenta of the 7r+ and the 7r~, which implies that the amplitude is zero. 
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represents the dominant contribution and the amphtude is proportional to — rud, which is usuaUy presented in 
one of the foUowing ratios 



that are connected via (r 



_ m,-m 1 _ ml -ml 

'7^£^ — rriu Q mi — rn^ 



f = ^R{r + l). (4) 



Consequently, a measurement of the decay rates of the — > Stt processes enables us a direct access to this difference 
(and by the use of the values m and from the lattice also to the individual masses of these two lightest quarks). 
Of course, in order for this extraction to be possible, it is necessary to have ChPT predictions for these decay rates 
with a sufficient degree of accuracy. 

Achie ving this is, however, a non-trivial task. The tree-level predictions, which are equivalent to the PCAC results 
(e.g. |16l4l8l |). would indicate a very large difference between m„ and md- Furthermore, the true energy dependence of 
the amplitudes is definitely different from the trivial one that PCAC proposes. The sizable one-loop corrections [l^ 
were still not sufficient to correct these discrepancies. At last, the inclusion of the two-loop corrections [l^l, which 
are also sizable, led to the predictions for both R and the Dalitz parameters describing the energy dependence of the 
amplitudes (cf. Tables HI and HIl below) that were in reasonable agreement with expectations. 

Nevertheless, if we study these results in greater detail, we find some hints that the 2-loop ChPT determination of 
R, which we are interested in, can still be inaccurate. The feature often put forward in this respect is the discrepancy 
between the experimentally measured and the predicted values of the Dalitz parameters (defined in Sec. Ill BT) . mainly 
for the neutral parameter a of (|27p . For a better quantification of the difference between experiment and theory, let 
us introduce 

2 _ / exp- theory \ 2 
^ = ( .[exp] ) ' 

where the theory enters only via its central value. Using this quantity when comparing the prediction of [20|. a — 0.013, 
with the best measurement of this observable by MAMTC j2i|, a = — 0.032±0.003, we obtain indeed a huge difference 
of ~ 225. However, there is a parameter for which this discrepancy is even more apparent, namely b from (|23p . 
Comparing the ChPT value b = 0.394 with the measurement by KLOE [l!] 6 = 0.124 ± 0.012 produces ^ 500. 
This raises the question about the origin of these discrepancies, and whether and to which extent they can also affect 
the determination of R. 

As was already stressed in (20| , the explanation of this difference between theory and experiment can be provided 
by the large theoretical error bars presented there (thereby making the theoretical and the experimental values 
compatible). The non-renormalizability of ChPT represents a major drawback of this theoretical framework when it 
becomes necessary to include higher and higher orders. Indeed, including two-loop effects to the computation means 
a rapid increas^ of the number of a priori unknown low energy constants (LECs) that have to be estimated before we 
can get a reliable prediction. We are far from a determination of all required LECs from experiment (or lattice) , and 
hence for many of them we have to rely on some estimates, predominantly of the resonance saturation type [23- 26,] . 
This brings an unknown error into the game — the error presented in [13] is an estimate by the authors obtained by 
taking the uncertainty of the amplitudes equal to one half of the two-loop contributions. 

Both the Dalitz plot parameter discrepancy and this drawback of ChPT affecting the predictivity of the chiral 
computation contributed to the development of alternative app roaches, among others the dispersive methods (27l - [29| 
and the non-relativistic effective field theory (NREFT) |30l434j . In order to understand their relative advantages and 
disadvantages, let us recall a few basic properties they share. All these approaches are constructed as effective field 
theories that on the basis of some assumptions (usually represented by some expansion of the amplitude) divide the 
phase-space of each amplitude into the "low-energy part" that is included in the computation and the "high-energy 
part" that is not known or at least less known. At tree-level one simply uses the amplitudes only in the low-energy 
region and is not concerned by what lies above the cutoff. In order to work consistently one needs to introduce 
a mechanism that picks up the contributions that contribute with the same importance to a given order, usually 
represented by a power-counting. Then, when computing the amplitude to the higher order, one needs to include also 



^Note that 0{p^) ChPT contains 102 (2+10+90) free LECs. In order to make this theory at the given order predictive, we would thus 
need to make at least 102 measurements determining these constants. Obviously, not all of these constants appear in a given amplitude 
— only a subset of them contributes to — > Stt, see below (Sec. Illlll . 
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loop contributions (either by means of taking into account loop Feynman diagrams, as a unitarity contribution, or by 
any other method), where one has to integrate also over the high-energy part of the intermediate amplitudes (over 
higher momenta of the intermediate virtual particles) . By using the "power-counting mechanism" or by adding some 
further assumptions, part of these contributions are considered negligible, but there always remains a part that is 
finite and unknown and has to be parametrized somehow — usually there occur new effective parameters in the model 
and the old ones are renormalized or shifted. Note that in ChPT that represents a Lagrangian effective field theory 
the power-counting mechanism is given by the chiral counting, which also monitors the number of LECs (effectively 
containing the contribution of the physics above the chiral cutoff — the hadronic scale) appearing at a given order. 

The importance of the one-loop (and in the recent years, also of the two-loop) tttt rescattering corrections has 
led |27h29|| to abandon, in a certain sense, strict chiral counting, instead attempting to obtain the amplitudes with 
two-pion rescattering effects formally included to all orders. These approaches employ a restricted version of unitarity 
(taking into account just the two-pion intermediate state), in the context of dispersion relations, the aim being to find 
a numerical fixed point solution of them. The mechanism assigning the importance to a given contribution is therefore 
based on the assumption that the two-pion rescattering effects are dominant. In the low-energy part of the amplitudes, 
the unitarity contribution of the physics above the threshold, where further intermediate states contribute and where 
the S and P partial waves of the considered amplitudes cease to be the dominant ones, are taken into account through 
subtraction constants. However, in order to restrict their number to a reasonable amount, one needs to impose some 
assumptions on the high-energy region (of both the physical amplitudes and of the amplitude constructed iteratively 
by the numerical method). In 28, 29] these assumptions are specified by the requirement to have only fouiH of them. 

The methods based on the modified non-relativistic effective field theory (NREFT) (30l - [3^ implement instead 
of the usual chiral expansion a combined expansion in powers of a formal non-relativistic parameter e and of a 
formal partial-wave tttt scattering-characteristics parameter a (representing scattering lengths and higher threshold 
shape parameters). The amplitude is then computed to the two- loop order in the NREFT Lagrangian formalism. 
The power-counting scheme is therefore based on the non-relativistic expansion together with the loop expansion 
(equivalent in this case to the expansion in the pion scattering parameters). In [s^j the results are presented including 
the orders up to 0(6*), 0{ae^), 0{a^e^), and partially also 0{a^e^) and 0{a^e^). By assuming that the included orders 
are dominant, the contribution of the intermediate states other than the two-pion ones have to be included through 
fouJ^ parameters coming from local interaction terms. 

Naturally, the reasonable question that has to be addressed in the future is whether each set of assumptions (either 
of ChPT, of the dispersive approaches, or of the modified NREFT) adequately describes the physics, and whether 
a possible drawback in this respect in any of them is paid off by the other advantages it possesses. The advantages 
and the disadvantages of these approaches were nicely summarized in [3^ . We emphasize only that NREFT provides 
analytic results that are easy to extent beyond the m^± = m^o limit, while the dispersive methods proceed numerically 
and their extension to full isospin breaking was never studied. On the other hand, whereas the NREFT expansion in 
e is safe only inside the Dalitz region, the results of the dispersive approaches should work also in some larger regions 
beyond it. Both of them have in common that, in contrast to ChPT, they directly use the physically measured tttt 
scattering parameters as inputs, but there remain fouJ^ free parameters that have to be fixed either from matching 
to ChPT or to experimental data. Moreover, these decays depend on i? or Q merely just through the normalization, 
which is factored out in both of methods. Thus, even if those representations are fitted to experimental data, the 
determination of i? or Q would still require to match with ChPT at least at one point, thereby fixing the normalization. 

The matching is not an easy task in this context. In addition to the differences in the structures of these results, 
since we are matching two different approaches with different power-counting schemes and assumptions, we need to 
find the region (or as discussed above at least one point) and the appropriate orders in both approaches in which 
their results are compatible. Nevertheless, thanks to the easy form of the one-loop ChPT amplitude, and to the fact 
that the physical reg ions of rj — )■ Stt decays are quite small, in both approaches the matching to one- loop ChPT was 
obtained (cf. [H, [23 and [3^). 

In conclusion, in order to determine the correct value of R from the r/ — !■ Sir decays, one cannot avoid discussing 
either the accuracy of the ChPT result for the amplitude and its possible corrections (by correcting the values of the 
0{p^) LECs Ci or by inclusion of some higher-order corrections into the ChPT calculation), or the existence of at 
least one point (or some region) where the current chiral result reproduces well the complete physical amplitude. For 
instance, the discussion of the influence of the C^s on the results can be studied using directly the ChPT amplitude, 
but its complexity and its extreme length together with the fact that it includes many such C^s complicates such an 
analysis. 

In [36| we have worked out a method using the dispersive relations and perturbative unitarity (i.e. a dispersive 
approach) for the construction of a representation of the 77 — )■ Stt decay amplitudes. We have given up on including 



Here we classify the number of parameters appearing in the case we take masses of the charged and the neutral pions equal, m^± = m^o- 
Note that from Sec.|V]it is obvious that in two-loop ChPT results there occur at least six independent combinations of LECs. 
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the TTTT rescattering contributions to all orders, but have instead required to obtain an analytic representation and 
paid care to the assumptions we are using in the construction, thus ensuring that the ChPT result can be obtained as 
a special case of our result. The method is based on very general principles, unitarity, analyticity, crossing symmetry, 
and relativistic invariance, combined with chiral counting. The fact that we require a representation valid to two 
loops in the chiral counting picks up the contributions that have to be included into the computation and tells us that 
at the low-energy region up to this chiral order all the other effects are taken into account effectively in terms of six 
subtraction parametersP. 

The full strength (and our original motivation) of this method arises when we want to include the isospin breaking 
induced by the mass differences between mesons belonging to the same isomultiplet (cf. 36]). However, even in the 
case where we consider the leading order in the isospin breaking, for which the two-loop ChPT result is available, 
the representation constructed by this method can be useful. Thanks to its simple and compact analytic form, to 
its capabilities to include all the chiral 0{p^) effects important in the kinematic decay region of 77 — > Svr into those 
six subtraction parameters, and to its easy correspondence to the ChPT, this representation is helpful when one is 
addressing the questions we have premised above, namely, whether one can obtain a reasonable agreement in the 
determination of the Dalitz parameters from experiment and from the NNLO ChPT amplitude with the corrected 
set of the C^s; how such a change would influence the determination of R\ possibly also whether there exists another 
simple way how to solve that disagreement. 

In addition, we do not need to work only in such a close connection to the two-loop ChPT amplitude. Our 
representation is more general than the two-loop ChPT amplitude (simply stated in the way that the values of our 
parameters need not to be held at the values stemming from the ChPT), based only on the specific chiral orders 
of the partial waves of the amplitudes (cf. e.g. (131). In order to respect such chiral power-counting, we need to 
distinguish between various orders of our parameters. By weakening this requirement and by a simple change of their 
interpretation we can perform a partial resummation that mimics a part of the previous dispersive approaches. By 
that we have therefore replaced the assumptions represented by the chiral counting with the assumption that the 
contributions we have included by this resummation are the dominant one. 

In any case, such representation is suitable for fitting the experimental data. We can thus change completely 
the strategy and instead of trying to correct the amplitude stemming from ChPT, we use our representation as a 
parametrization of the data, from which we can compute the value of R. However, as was discussed above also in this 
case, we need to fix the normalization from ChPT. For that end we need to find a region where the chiral expansion 
of the amplitude converges fast, where the two-loop ChPT amplitude reproduces the physics well. Thanks to the 
form of our representation and its simple connection to ChPT the analytic dispersive method is helpful also in this 
analysis, resulting with the recipe for such matching. 

We want to stress that in [36j the inclusion of the isospin-breaking corrections stemming from m^± ^ m^o is 
presented, but as we discuss in Sec. I VIII the current experimental data do not yet allow to perform the isospin- 
breaking analysis. Thanks to the planned improvement in the neutral decay measurements (cf. [s^) we should add 
that this possibility is just behind the corner. In the paper we therefore work in the limit m^ri = TOtt" with the 
exception of a few discussions of the effects appearing beyond this limit. This discussion in full detail is however 
planned in our next paper [soj . 

Through relation ([20]) this limit connects the charged rj — Tr+Tr"??" decay with the neutral one, rj — 37r°. Using 
our representation on the basis of the above mentioned analyses of the charged KLOE data [2^ (the most precise 
measurement of this process that exists) , we can therefore determine the values of the neutral Dalitz decay parameter 
a and discuss its connection to the direct neutral measurements (from Table lll| . 

The plan of our paper is as follows. After recalling some basic properties of the amplitudes of these decays and 
introducing our notation in Sec. [Hi we recall in Sec. Illll the ChPT computation of 77 — ^ Stt with the special emphasis 
on the contribution of 0{p^) LECs to the Dalitz plot parameters. From that analysis there follow a few combinations 
of observables that should be (at least in the first approximation) safe from the incorrect determination of these LECs. 
In Sec. IIVI we present the dispersive construction of our representation for the 77 — >■ Stt decay amplitudes. Section IVl 
discusses the connection between our representation and the ChPT result. Section IVIl is then devoted to the numerical 
analysis of the charged decay. We start with the determination of the values of our parametrization for NNLO ChPT. 
Then, inspired by the result of Sec. IIIIl we study the influence of changing the C^s in the NNLO ChPT amplitude in 
order to reproduce the charged KLOE data !22] on the values of the physical observables we are interested in. Then 
in Sec. I VI CI we perform the analysis in which the values of the dispersive parameters are set by KLOE and only the 
normalization is determined from ChPT. In that section we present also the procedure of the matching that should 
reduce the uncertainty coming from the chiral expansion of the amplitudes. In Sec. IVIII we use these analyses also 
for the determination of the neutral Dalitz decay parameter a and discuss briefly the determination of the ratio of 
the neutral and the charged decay width. Finally, our conclusions can be found in Sec. I Villi We have devoted two 
appendices to the discussion of properties of the kinematic functions appearing in our dispersive representation. 
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II. BASIC PROPERTIES 
A. Kinematics and notation 

We are interested in two decay channels of rj, the charged one rj — tt+tt^tt'' and the neutral one t] — > 37r°, generically 
denoted as 

r^{P)^7:^ip,)7:\p,)7r°{p,). (6) 

The amplitudes of these decay processes can be obtained by analytic continuation of the amplitudes of the corre- 
sponding scattering process 

viP)n°ip) ^ n^{pi)7r\p2) (7) 

by taking p = —ps. In the following sections, we use the usual Mandelstam variables. In the scattering region they 
are defined by 

s^(P+p)\ t^{P-p^)\ U^{P-P2)\ (8) 
while in the decay region we take 



{P-PjY- (9) 



These variables satisfy 



where s = t = u = s^, with 



t + W = Ss'^ = Sl +S2 +S3, (10) 



35^= = to; + 2m| 2 + (11) 

corresponds the center of the Dalitz plot. Here TO3 = m^o, whereas toi^2 = itl^q for rj — > Stt", while mi^2 — '7ijr± 
for 77 — > vr+Tr^Tr". Up to a convention-dependent phase factor, the crossing relation then means a substitution of the 
variables (s, t, u) by (53, si, S2), together with the appropriate analytic continuation from the scattering to the decay 
region. Bearing this in mind, we can therefore interchange freely between these two sets of variables. 

The constraints (|10l) tell us that just two of the kinematic variables are independent. We can choose them to be, 
for instance, S3 = s and si — t. The plot of the dependence of the decay amplitudes on these variables is called Dalitz 
plot. The physically allowed kinematical regions for the different crossed amplitudes are constrained by kinematical 
limits arising from the condition that the energy of a real particle has to be at least equal to its rest mass. Therefore, 
for a decay process the variable S3 is bounded by 

(toi + TO2)2 < {pi + P2f - S3 ^{k- psf < (to^ - m3)2, (12) 

whereas for a scattering in the s-channel 

s > (m,, + TO3)^. (13) 

For a fixed value of S3, we obtain bounds for the physical values of si (and similarly for S2), {ss) < si < 3^(33) 
with 

± 3s'=-S3 A^3Ai2 , A;^f (S3)A}/^(S3) 

^''^ = —2 2^ ^ 2^^ ' 

where j = rj, 1, 2, 3) 

Ay- = to| - m], (15) 

Aij(s3) = (S3 - (TOi + TOj)^) (S3 - (jUi - TOj)^) . (16) 

Since for both cases under study m\ = m2, the bounds simplify to 

(35^-53 ±Aif(s3)a(s3)), (17) 
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with 

/ 4rn? n 

^(53) = Jl -■ (18) 

V S3 

As was recalled in the Introduction, the amplitudes of the processes 77 Stt, A{s,t,u), are proportional to the 
difference of and m^. We therefore pull out this factor, defining 



\/3 1 1 

Ais,t,u) = —-j^M{s,t,u), (19) 



where the ratio R, which is defined in ([3]), measures the relative violations of SU{3) and of SU{2), and F^^ is the 
physical pion decay rate (in our numerical analyses of Sec. IVII in order to be in correspondence with [2C|, we tak^l 
= 92.4 MeV). In accordance to the notation introduced in our general paper [s^, when the distinction becomes 
necessary, the quantities associated to the charged (ry tt+tt^vt") or neutral decay {rj — )■ Stt") are denoted with the 
subscript x or 0, respectively. 

In this paper we work mainly to lowest order in the isospin breaking, i.e. we consider the case where all isospin 
breaking is contained already in the normalization prefactor R from ([T5| , and the rest of the amplitude is computed 
in the isospin limit. Then due to the isospin structure, the amplitudes A4{s,t,u) are related by 

Mq{s, t, u) = -Mx{s, t, u) - Mx{t, u, s) - Mx{u, s, t) (20) 

(the minus sign is due to the Condon and Shortley phase convention) and in both A^o('S,i,w) and AAx{s,t,u) there 
appears only one pion mass m^. It is why we refer to this case as the m^± = m^o limit, or more loosely as the isospin 
limit. However, when making comparisons with the ChPT calculation of [20| . we use exactly the same values for rj 
and TT masses as were used there[£|, and take for the isospin mass in each process a different value — in the case of 
the charged rj decay we take 3mJ = 2m^± + m^o , whereas m^r = m^o in the case of the neutral decay. So defined 
has the advantage that in both processes we reproduce exactly the physical location of the center of the Dalitz plot 
and reproduce almost exactly the physical value of the normalization Q,^ of the Dalitz variables x and y from (j2ip 
below. When computing the integrations over the phase space used for setting the normalization from the measured 
decay rate, we employ again the physical m^ri and 711^^0 masses for the determination of the phase space. 



B. Dalitz plot parametrization 

The standard parametrization of a decay process P -> Btt is called a Dalitz plot parametrization (cf. [llj). It is a 
polynomial expansion of |-4(si, S2, 53)^ around the center of the Dalitz plot. The parameters are usually normalized 
in order to be dimensionless. The variables of standard use for the charged rj decay are then 

^^^73 ^1-^2 ^ \/3 (S2 - Si) y^^_2=— ^((771 - 77^ o)2 - S3) - 1 (21) 

where Tj is the kinetic energy of the j-th pion in the 77 rest-frame. For the charged decay the energy of the reaction 
Qri = 771,, — 27n^± — uIt^o whcrcas for the neutral one Qri — rriji — 3m^o. In the case we use in this definition the 
physical values of the masses, for the charged decay the point x — y = 0, around which we expand the amplitude, 
does not coincid^ exactly with the center of the Dalitz plot. However, in the isospin limit, 

and the center of the expansion x = y = matches the center of the Dalitz plot. 
The parameters relevant to the decay 77 — > 7r+7r~7r° are usually labeled according to 

\Mx{si,S2,S3)\^ = \A\^{l + ay + by^ + cx + dx^ + exy + fy^ +gx^y...) , (23) 



^Recent analyses (e.g. 140(1 ') indicate a slightly smaller physical value, = (92.22 ± 0.07) MeV. In order to fully include this change into 
our computation, redoing of the analysis 20j with new values for and for the pseudoscalar masses would be necessary. Note, however, 
that a mere change of this value just in this definition leads to a shift of the value of R of about 0.4%, which is negligible with respect to 
the other sources of errors occurring in the presented results. 

^Note that the value used rrirf = 547.3 MeV differs slightly from the current PDG value [4l| . 

^The X = y = point is slightly shifted in the S3 direction to S3 = s'^ + |(m^± — m,^o)(2m,, — m^± — m^o). 
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Gormley et al.[42j 


-1.17 ±0.02 


0.21 ±0.03 


0.06 ±0.04 




Layter et al. [43] 


-1.08 ±0.014 


0.034 ± 0.027 


0.046 ±0.031 




Crystal Barrel [44] 


-1.22 ±0.07 


0.22 ±0.11 


0.06 ±0.04 




KLOE [22] 


-1.090 ± 0.020 


0.124 ±0.012 


0.057 ±0.017 


0.14 ±0.02 


ChPT NNLO [20] 


-1.271 ±0.075 


0.394 ±0.102 


0.055 ±0.057 


0.025 ±0.160 



TABLE I. Values of the Dalitz plot parameters of the 77 — >■ tt^tt tt decay coming from various experimental and theoretical 
determinations . 



where A is the value of the amplitude M.^ at the point x = y = Q. Charge conjugation forbids the appearance of 
terms containing odd powers of x in this expansion, and so c = e = 0. 

The values of the parameters obtained by various experiments are listed in Table ID These values are compared with 
the NNLO calculation in ChPT f20|. All of the experiments find the values of c and e compatible with zero. From 
the table it is obvious that the precision of the determination from KLOE exceeds significantly the precision of 
all the other experiments, which are more than ten years older. It is also up to now the only experiment that has 
determined the parameter / with a reasonable precision. 

At this point let us also mention the linear Dalitz parametrization for the amplitude itself (cf. Appendix A of 20] ) : 

M^{s,t,u)= A{l+ay + by^ + dx^ + fy^ + gx'^y + ...), (24) 

where the parameters a,b, . . . can now be complex in general. (We have already omitted the terms violating the 
charge conjugation symmetry of the amplitude.) The parameters of (|23p can be expressed in terms of these linear 
Dalitz parameters — the relations are simple to obtain by squaring (|24p and by comparing the terms with the same 
powers of x and y. 

At leading order, the parametrization of the 77 — >■ Stt" differential decay rate depends only on the kinematical 
variable 



3(.g)2) , (25) 



which denotes the distance from the center of the Dalitz plot, normalized to one at the edge of the decay region. 
However, higher orders corrections do not preserve this accidental rotational symmetry, and we need again x and/or 
y from relations ((2T|) . Note that the relation 



z = ± y2 (26) 
holds. The Dalitz plot parametrization for this process reads 

, . = 1 ± 2az + 2/3y(3z - 4?; ) + 72 + . . . . (27) 

The factor of 2 in front of a and /3 is a mere convention to stress the connection with the direct linear Dalitz 
parametrization of the amplitude itself (see (|24p above). For a better visualization of the violation of the rotation 
symmetry in {x,y)- plane at higher orders, it is convenient to introduce the polar coordinates (cf . also [3^ ), x = p cos 0, 
y — psn\4> with distance — z, for which we have 2/(32: — 47/^) = p^ sin(3(/)). 

Various experimental and theoretical determinations of the parameter a are given in Table [iTl Note the sign 
discrepancy between the ChPT determination (with however large error bars) and the experimental measurements, 
which we will briefly address in Sec. I VIII Up to now, no experiment has so far published any constraint on the other 
parameters, such as (3. 

In the case we work to first order in isospin breaking, the isospin relation (|20p translates into the following relations 
between the neutral Dalitz parameters and the parameters of the linear parametrization (I24p (cf. again Appendix A 
of HI) 

a = i(Re5 + Re J), (28) 
/3=i(Re5-Re/). (29) 
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a 



Crystal Barrel [45] 


-0.052 


± 0.020 


Crystal Ball [46] 


-0.031 


± 0.004 


WAS A/CELSIUS [47] 


-0.026 


± 0.014 


WASA/COSY [48] 


-0.027 


± 0.009 


Crystal Ball @ MAMI-B [49] 


-0.032 


± 0.003 


Crystal Ball @ MAMI-C [21J 


-0.0322 ± 0.0025 


KLOE [50] 


-0.0301 ± 0.0050 


ChPT NNLO [20] 


0.013 


± 0.032 



TABLE II. Experimental and theoretical values of the slope parameter a of the 77 — s- 37r decay. 



They can be rewritten in terms of Dalitz parameters of the charged decay. However, there still remains a dependence 
on the imaginary parts of the linear parameters, 

\ ( „ a,, „ c?\ Ima ( - - Ima\ 

'g_/+_(6_rf)__j +_ flni6-Imd-a— j . (31) 



C. Adier zero 



The isospin-breaking part of the QCD Hamiltonian density ^ can be written as (in this subsection A"^ are Gell- 
Mann SU(3) matrices) 

^qcd(2;) = (m-d - mu)Sz(x), (32) 

where 

S^{x) = \^{x)\\{x). (33) 

Therefore, to first order in rud — mu, the ampUtudes of the isospin-breaking processes that are described by this 
Hamiltonian are connected to Green functions with one insertion of zero momentum scalar density S3 (calculated in 
the limit ttiu = nid)- In the SU{2) x SU{2) chiral limit m„ = ma = 0, pions are genuine Goldstone bosons. For the 
corresponding amphtudes (/ -I- 7r"(p)]S'3(0))i) with a pion in the final state, we can thus derive the SU{2) x SU{2) 
soft-pion theorem in the general form 

(/ + ^"(p)153(0)li) ^0 for p->0. (34) 

The derivation of the theorem proceeds in the usual way, except that now, because of the insertion of 6*3(0) trans- 
forming under the axial SU (2) rotation nontrivially as 

5tS\Q)^'-q{x){\\\^}^^q{x), (35) 

it only holds provided a — 1,2. (For a = 3 there occurs an additional contribution from 5|5''^(0), which does not 
vanish.) 

For the decay 77 — > 7r+(p+)7r^(p_)7r°(po), this means that the amplitude Mx{s,t,u) defined in (|19p vanishes in the 
SU{2) X SU{2) chiral limit for either — Q or p- = 0, i.e. it develops two Adler zeroes [HI, [13] s = u = t — 
and s = t = u = rn^. As a consequence, expanding the amplitude Mx{s^t,u) beyond the SU{2) x SU{2) chiral 
limit in the independent variables s and (t — u) around the points where s = 0, (t — uY = m^, or more specifically 
around the points 

s = 0, 

^= ^("i^ + 3m^±m^), (-3g>) 
w= ^(™^ + 3m^=Fm2), 



according to (here we use the t -(r^ u symmetry of the amphtude) 



M,{s,t,u)^ J2 c^Js\{t~uf ~mt,y, (37) 

ij>0 

we can restate the above SU{2) x SU{2) theorem in the form 

hm Coo = 0. (38) 

niu ,md— >0 

Since the position of the Adler zero is determined up to 0{rn^) corrections, an analogous statement remains true also 
for similar expansion coefficient Cqq''''' corresponding to an expansion around the points with s ~ S.m'^, (t — u)^ = 
(m^ + Cm^)^, namely, around the points 



s 



H^X) ^t^2^CT0ml, (39) 
"(«,C) ^'"T i(C±0™^> 

where ^, C are reasonably small and behave as 0(1) for mu,md — ?> 0. For the value of the amplitude at these points 
we therefore obtain 

■^-(s(«.C)'*(«,C)'"«,C)) = 4o'^^ - 0(774) (40) 

and its absolute value is expected to be numerically small. 

Note, however, that the remaining coefficients Cij are not protected by such a factor m^, and the same is true also 
for the value of the amplitude at points far from (s, t, u), where M^is, i, u) can be enhanced by a factor m^/mj with 

respect to A4x{s,t,u). Note also that a small numerical value of cqo (or Cqq'^"* in general) does not necessarily imply 
that its chiral expansion shows better convergence than the one of any other Cy , in the sense that for the ratio of two 
subsequent corrections the relation 

(„+l) (n+l) 
t-nn 



« ^ (41) 



■-00 

does not necessarily hold. 



D. Isospin violation and cusp 



In the case we go beyond the first order in the isospin breaking, in addition to the complications that the two 
?7 — > 37r decay amplitudes are no longer connected by (j20p . and that the expressions are more complicated due to the 
fact that there appear two different masses of pions, in the processes with two neutral pions in the final state there 
occurs an interesting phenomenon called cusp. This effect is caused by different charged and neutral pion masses and 
is connected with the contributions of tt+tt" intermediate states rescattering back to 7r°7r°. Such a state generates a 
square root singularity, which resides at 4m^±, lying above the physical threshold, 4m^o, and the unitarity cusp is a 
result of the interference between the part of the amplitude containing this singularity and the rest without it. 

It is obvious that the cusp emerges only in the case when isospin breaking is included also in A4{s, t, u) and that its 
strength is sensitive to tt+tt" — )• -k^tt^ scattering at threshold (mainly to the scattering length of this process). This 
property can be used for a determination of the scattering lengths from the measurement of the cusp |53l - [55| . 

Let us try to estimate the relative sizes of the cusps in various processes where a pseudoscalar, namely Kl or 
ry, decays into three pions. (This discussion is inspired by [5g] and [57|.) Because the pion rescattering part will be 
approximately the same for all the processes, we may consider the notion of "visibility" of the cusp in these processes 
by comparing the relative ratios between the cusps and the regular parts of the amplitudes. 
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where \A\ is the absolute value of the matrix element of the indicated process and Kc is a multiplicity factor corre- 
sponding to that process, equal to 2 in the case the decaying particle is charged (two possible tt+tt" — > tt^tt^ scatterings 
are then possible), and to 1 in the other cases. These ratios have to be evaluated at the cusp point s = 4m^±. 

Using the measured relative decay rates and the values of Dalitz parameters from [4l|, we obtain for these processes, 

7(i^±)-7.3, 7(ifi)-0.45, 7(7,) ^ 0.93. (43) 

From that we can estimate that the effect is approximately 16 (8) times more pronounced in the decay with 
respect to [rf) decay. 

First indications of the cusp effect were already observed also in the -q — >■ 37r° decay (cf. e.g. [21]). This effect 
however appears at the edge of the decay regiorQ and is therefore not simple to measure. 

For the time being, because of this lack of data, we shall not pursue the discussion about the cusp here (even 
though our representation describing also this effect is prepared [36|, HI]). Instead we will work in the strict isospin 
limit beyond the trivial order at which 77 — !■ 37r decay is forbidden, i.e. AAx or AAq is taken in the isospin symmetry 
limit. 



III. CHIRAL PERTURBATION THEORY 



Let us briefly recapitulate the ChPT calculation of 77 — > 37r decays. As was discussed in the Introduction, direct 
electromagnetic corrections to these processes are very small, and thus they proceed mainly through the part ^ of 
the QCD Lagrangian. The leading order (LO) calculation was performed in [16l - ll8l |. which in our notatioro takes the 
very simple form 

M^^\s,t,u)^^ml-s. (44) 

The next-to-leading order (NLO) was provided in |T^. Its form is discussed in Sec. lIV Al below. The 0{p^) corrections 
were studied quite recently in . From these three successive orders one can see that ry — ^ 37r thus represents a case 
where the chiral corrections are larg e 1591 . Moreover, it seems that also the two- loop ChPT result supplemented with 
the existing LECs determination of [20l| is not working very well as we have demonstrated on the example of Dalitz 
parameters in the Introduction. 



A. Contribution of the constants d to Dalitz parameters 

In the NNLO result there occurs a great deal of 0{p^) low-energy constants C,; which are only estimated from 
resonance saturation. Hence, the first question one has to ask is whether the discrepancy with experiment cannot be 
accounted for by the unsatisfactory knowledge of the 0{p^) low-energy constants. 

Let us thus study the contribution of 0{p^) LECs (Cf ) to the Dalitz parameters of the individual decay modes. 
There are several possibilities how to determine these parameters from the computed amplitude Ai{s,t,u). For 
instance, we can expand \A4{s,t,u)\'^ to the order 0{p^), and then make the Taylor expansion at the center of the 
Dalitz plot. This would result in the linear dependence of the Dalitz parameters on the C^s. Provided we did not 
chiral expand first and instead made a fit of the modulus squared of the complete 0{p^) amplitude to the Dalitz 
parametrization (as it was done in i2C|), we would get a more complicated dependence on the C^s including also 
quadratic and mixed terms. Such contributions should be, however, suppressed by the chiral counting. Nevertheless, 
they can bring sizable changes in the final numerical predictions. In order to obtain the linear contribution only, we 
follow the first procedure. 

We start with the neutral decay mode. The explicit dependence of a on the C^s was already given in [33 |. 

16m2(77i,, - 3m^)2 
ac = (45) 



^In the {x,y) plane, the cusp is located on the segment y = yiim? ^) ~ 0.773 and on two other segments obtained by s •<-)■ t and s -o n 
(i.e. obtained by rotation of the original one by ±120° around the center of the Dalitz plot). Its position thus does not respect the 
accidental rotation symmetry, and depending on its direction in the {x,y) plane, the corresponding value of z changes from 0.597 to 0.883 
as — 7#^, with 00 = 0°,±120°. 

"In this work for the various chiral orders we follow the convention of 12011 . where the amplitudes are at a given order simplified using 
Gell-Mann-Okubo relations, physical decay constants and physical pseudoscalar masses and the so induced differences are included into 
higher orders. 
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with 

C„ = a, + CI + iCl + C[o - 2Cl^ + 2Cl^ + 36*2^4 + C^^^ . (46) 

Further, by a careful investigation of the 0{p^) polynomial of the amplitude calculated in [2^, we realize that there 
is no contribution of the ds to the second neutral Dalitz parameter (it is connected with relation (|29|) and the fact 
that in the charged decay fc = gc as stated below) , 

Pc = 0. (47) 

In the case of the charged decay we summarize first the contributions of the C^s to the linear coefficients a, 6, d, 
/, g defined in (|24p that are directly connected with the amplitude. These parameters can be in general complex but 
since we deal only with the linear contribution of the C^s, they contribute only to their real values. 

By a simple algebra one obtains 

«c = - (KC - em^m^Cr - SmtC:) , (48) 

where we have slightly more complex structure 



C2 = C[ - 2C; - 6Cl - 2Cl - ZCl + + C[o + 6C[i - 20C[2 - ISC^g + 9C[4 - 36C[e + 9C[7 + 45C[8 

- 81C[9 - 54C2^o + 8C2^2 + 9C2^4 + ^2^5 - 18^^27 " 36^2^8 - 54CJi - 54C7J2 " lOSCg^g- 

C;"^ = C[ - 2Cl - 2Cl + CI + 2Cl ~ C[o - 4C[i - 4C[2 + 12C[3 - 4C[4 - iCl^ - 12C[e + 2C[7 + 6C[8 
+ 4C22 + 3C24 — C25 — 2C2g — 6C27 ~ 126*28 + 4C29, 

= C[ - 2C3" + 2Q' + 3Q + 3Q - 3C[o - + 12C[2 + ISCrg + 3C[4 - 6C[5 + 126^6 - 9Cr7 + 3C[8 

— 27C{g — I8C20 ~ 3C24 — 36*25 ~ 4C26 + 6C27 + 126*28 + 8629 — 186*34 — I8632 — 36633. 

Similarly, we have 

,2 



(49) 
(50) 
(51) 



8m,,(m„ — 3m7r) , ^ „ n \ , . 



where 



= 6[ - 26^ + 26^ + 4C; + 36^' + 68^ + 66^' + 6[o - 66[i - 26[2 + 186^3 + 26^2 + 66^4 + 46^5, (53) 
= 36[ - 663^ + 6Cl - 2Cl + ZCl + CI - 66^ - 56[o - 66[i + 106[2 + 186[3 - 106^2 " 66^4 + 46^5, (54) 

= - OF4r^2_^2^ y-^r^^d + ^^.^d ) , (55) 



3F4(m2 



with 



Cl = 6[ - 263^ + 264^ + 36^' - 368^ - 669^ - 36[o - 66[i + 66[2 + 186[3 - 66^2 - 662^4, (56) 
Q = 36[ - 66^' + 664" + 265" + 366^ + 56^' + 669" - 6[o - 66[i + 26[2 + 186[3 - 262^2 + 662''4 + 862^5, (57) 



and 



with 



16m^(TO^ — SniTr)^ 

■^c ^ qf4^ 2 2\~*^/ ^8) 

3F4(m^ - mi) 



Cf = C[- 2C; - 261 . (59) 
Finally, the contribution of 6^8 to the parameter g is the same as in the case of the parameter /, 

gc = fc- (60) 



12 





KLOE 


ChPT 


ChPTg 


NREFTi 


NREFT 


a 


-1.09 ±0.02 


-1.271 ±0.075 


-1.262 ±0.079 


-1.279 ±0.012 


-1.213 ±0.014 


b 


0.124 ±0.012 


0.394 ±0.102 


0.407 ±0.108 


0.361 ±0.021 


0.308 ± 0.023 


d 


0.057 ±0.017 


0.055 ±0.057 


0.081 ±0.089 


0.053 ± 0.003 


0.050 ± 0.003 


f 


0.14 ±0.02 


0.025 ±0.160 


0.009 ±0.170 


0.089 ±0.018 


0.083 ±0.019 


9 


~ 





-0.07 ±0.19 


-0.043 ± 0.002 


-0.039 ± 0.002 


a 


-0.030 ±0.005 


0.013 ±0.016 


-0.024 ± 0.005 


-0.025 ± 0.005 


P 




-0.002 ±0.025 


-0.004 ± 0.001 


-0.004 ± 0.001 



TABLE III. Dalitz parameters for charged and neutral decay mode of ^ 37r calculated in ChPT at order 0(p®) 0|, NREFT 
[3^1 (in addition to their final value given in the ultimate column, in the penultimate column captioned with NREFTi we list 
also the values without isospin-breaking effects included) and as measured by KLOE [13, H^. ChPTg represents a fit redone 
with g included (but without any higher parameters), cf. fit NNLOq in [20| . 

Using these relations in the same spirit as in (60j , we can construct the combinations of physical (or quasi-physical) 
quantities which do not depend on any Cc 

1. (6± J- 2a)|^ = ^ 

reli EE (4(6±d) - - 16a)|^ = (61) 

2. (/-5)lc = 0^ 

rel2 = (a^ -4a6±4ad + 8/-85)|p = (62) 

3. /3|c = 

4. Cl + Cl-C^,-Cl = Q 

Let us discuss them in more detail (in a reverse order). The last expression, of course, does not represent any 
combination of physical quantities, and so it is not possible to use it directly in connection with any observable. It 
could be, however, useful for lattice simulations, where one can vary the meson masses. On the contrary the third 
relation, stating that the second neutral Dalitz parameter j3 does not depend on any Ci, represents a simple possibility, 
open to future experiments, how to check the ChPT results unaffected by the error stemming from the estimates of 
Ci. Now let us turn our attention to the relations reli and rel2. The latter was mentioned in [60], while the first one 
was implicitly stated in In fact, reli is a simple consequence of the isospin relation ([50]) stating that the C^s do 
not contribute to Ima and should thus be valid for any real contributions to the Dalitz parameters appearing there 
(not only for the contributions of the Cis). 

For the comparison of ChPT results [lol with the values measured by KLOE [S^ we can use not only the values 
of the Dalitz parameters summarized in Table IIIII but also the combinations of these parameters (pT|) and ([5^ 
that are (at least in a first approximation) Ci-independent. It means that the influence of all physics beyond the 
pseudoscalar domain (mainly scalar and vector resonances) on these combinations is hidden in O(p^) LECs L^, which 
are phenomenologically much better under control than the C^s, thereby providing a clearer theoretical output. We 
should remark, however, that the independence of all these relations on the C^s occurs only in the case we take 
m^± = TO^o. Away from this limit the situation can be different and these combinations can still have non-negligible 
dependence on the 0{p^) LECs. 

The values of these combinations that use the data from Table IIIII are presented in Table IIVI This table indicates 
that even though the central values of the individual Dalitz parameters determined by ChPT and KLOE differ, the 
central values of these two combinations are in a good agreement, which indicates that ChPT is not working at 
all that badly. Unfortunately the large errors quoted there somehow put down the importance of any conclusions. 
However, one should bear in mind that the values quoted in Table IIVI were computed just using the values and 
the error bars presented in Table IIIII that were attributed mainly from the fitting procedures and are thus strongly 
correlated. This can affect the positions of the central values by small changes, but primarily the error bars of these 
combinations are then overestimated. Note that the errors of the Dalitz parameters from ChPT are enhanced also 
by large systematic uncertainties of the amplitudes entering these fitting procedures. Such uncertainties were caused 
mainly by uncertainties of the C^s, which should be substantially eliminated in these combinations. We also observe 
another artifact of the fitting procedure when comparing the values denoted by ChPT and ChPTg that differ just by 
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KLOE 


ChPT 


ChPTg 


NREFTi 


NREFT 


reli 


0.02 ± 0.12 


-0.03 ±0.72 


0.15 ±0.79 


0.41 ±0.12 


0.35 ±0.13 






0.1 


1.2 


10 


7.6 


rel2 


0.12 ±0.21 


-0.13 ± 1.4 


0.28 ±2.1 


0.54 ±0.19 


0.44 ± 0.20 






1.4 


0.6 


4.0 


2.4 



TABLE IV. Two relations defined in (|6ip and ((62]) compared using different models. The quantity defined in ((5]) was 
computed for the theoretical frameworks by taking for the experiment the values from exp = KLOE as given in the first 
column. 



the truncation of the Dalitz parametrization at / and g, respectively. The combination reli is according to relation 
(pO]) equal to 4(Ima)^, which should be therefore positive. The value denoted by ChPT does not possess this property 
even though that the value of 4(Inia)^ obtained by a direct fit of the original amplitude in [2^ reproduces well the 
value given in the column ChPTg. 

A similar effect can occur also for the KLOE values since in 22 1 the valu^ of g was not presented (only its 



compatibility with zero). As an illustration, we remind the reader that if we added to the values of a, b, d, f measured 
by KLOE the value of g = —0.02 {g — —0.04), we would obtain an exact match of the so defined experimental value 
of rel2 with the value from ChPTg (NREFT). 

In these two tables we have also studied predictions of NREFT [s^ . Since that method is built in a different way 
than ChPT, the combinations of the observables appearing in reli and rel2 have no special significance there. However, 
they are still valid combinations of observables and so nothing prevents us from using them for comparison of the 
predictions from any theory with the experiment. The lesser agreement of NREFT and KLOE in reli was already 
pointed out in [33| in terms of different values of Im a stemming from the representation of [s^l and the one coming 
from the KLOE measurement and the relation ([50]). Together with the slight inconsistency also in rel2 depending 
only on the parameters of the charged decay, this indicates that there is a problem either on the side of the current 
determination from the KLOE group or on the side of the NREFT representation. 

We conclude this discussion with the statement that a new measurement of the charged Dalitz parameters (possibly 
taking into account these two relations) would therefore be highly desirable. Before that, we are not able to answer 
the question whether it is possible to reproduce the physical Dalitz plot distribution with a better determination of the 
LECs Ci or whether the discrepancy between the ChPT-computed and the experimentally measured distributions has 
some other origin (slow convergence of the chiral counting, . . . ). In addition, should the experimental value confirm 
the values inconsistent with the pre dictions of ^34i] . even if one accepts the explanation for the discrepancy of the 
neutral parameter a proposed in [3J| , the issue of the discrepancy for the charged parameter b would remain open. 

But for now, inspired just by the quite good consonance of the current KLOE and the ChPT values of the Ci- 
independent relations, we would expect that by finding the right values of the CjS we would reproduce (at least 
partially) better the physical values of the Dalitz parameters. The natural question can arise now whether it would 
be possible to find an elaborate determination of such C^s going beyond the crude estimate of the simple resonance 
saturation model used in [23| • 

Let us start with a. Its resonance saturation is simpler as there are no vector resonance contributions. For the 
simple scalar resonance model used in [23| we obtain (cf. (sj) 

16rrii(m„ — 3'mT^)^CdCm 



that is positive. However, the minimal chiral symmetry breaking introduced in 61| changes Cm into €,„ + c^e 
and especially for standard hierarchy (e„i < 0) one can thus produce a negative contribution to a. Using the same 
numbers as obtained from the phenomenological study in (61j . where they distinguish two models, one representing 
the inverted hierarchy (the model called A) and one representing the standard hierarchy (called B), we obtain 

a\c [A] - 0.002, or a\c [b] = -0.0005 , (64) 

which lead to the final values a = 0.01 and a — 0.007, respectively. 

The situation for the charged decay mode looks more complicated. The transition from the amplitude to the Dalitz 
parametrization leads to many mixing terms and the dependence on the C^s is not linear. Even though, as already 



^Note the different notation of this Dalitz parameter in [22| — for the parameter denoted in this text by g, KLOE uses symbol h. 
"Naturally, repeating the KLOE fit with g included would also change the values of the further parameters (cf. again the difference between 
the values from ChPT and ChPTj,). 
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model 


\Af 


o 


b 


d 


/ 


9 


ChPTg 


534 


-1.26 


0.41 


0.081 


0.009 


-0.072 


simple 


516 


-1.39 


0.47 


0.10 


0.025 


-0.088 


model A 


723 


-1.31 


0.41 


0.081 


0.024 


-0.069 


model B 


1835 


-1.19 


0.33 


0.052 


0.020 


-0.040 



TABLE V. The comparison of the impact of the different models for resonance saturation on the Dalitz parameters. For vector 
resonances we have taken the model from [6^ . different cases for scalar resonances are: the simple one from [20| ]. and the models 
A and B from 

mentioned, such higher terms are theoretically suppressed by chiral counting, in practice they can turn out to be more 
substantial than anticipated (it is true especially for model B). In order to get more reliable results we perform a full 
fit to the Dalitz distribution in exact correspondence with j^Oj , with the exception that we fit a polynomial of the third 
order (i.e. including g), which corresponds to ChPTg in Table llTll The vector resonance saturation employed here is 
based on the model and the phenomenology constraints from j62| . The resulting fits to the Dalitz parametrization 
are summarized in Table |Vl It is clear from this table that, as we have anticipated, the CiS have a bigger effect than 
expected from mere chiral counting. They also have an impact on the normalization which in the case of model 
B is far from being negligible. Let us note at this point a few things concerning the resonance saturation. It is 
obvious from Table IVl that model B would produce an unrealistic increase of the amplitude (thereby also of i? or Q). 
It does not, however, mean that this model for scalar resonances is disqualified. Higher resonances, representative of 
the physics beyond the pseudo-Goldstone bosons, contribute to both LjS and Qs (when talking about NNLO). One 
cannot just keep their influence on C^s ignoring their presence in L^s and thus merge inconsistently two models, i.e. in 
our case the model used in "flt 10" of [63] and the models A or B. One can always try to be as "harmless" as possible 
with any extension of the simple resonance saturation and try to preserve the original values of LjS (as was to some 
extent possible for the chiral symmetry breaking construction done in [S^l), hoping that the new effects induced by 
the new resonance terms will not change considerably the original and phenomcnologically successful "fit 10". But 
generally this is not guaranteed. 

The detailed analysis based on the current experimental data which would take into account simultaneously and 
consistently various resonance estimates for both O(p^) and 0{p^) LECs is beyond the scope of this paper (however 
such a project is under investigation [64','65']). Instead we present another representation that can be used for analyzing 
the data without addressing the values of the individual CiS. 



IV. DISPERSIVE CONSTRUCTION 

The dispersive construction to be presented below is based on the reconstruction theorem [66l - [68t , which takes into 
account only the most general properties of the amplitude, namely, relativistic invariance, unitarity, analyticity and 
crossing, supplied with chiral counting (e.g. expansion in powers of momenta and of masses of the pseudoscalars). 
This framework provides the most general form of the amplitude under consideration in the low-energy domain, up 
to a remainder of the chiral order 0{p^). Such a construction requires at the same time the scattering amplitudes 
related to the original one by two-particle unitarity and by crossing. (Contributions to the unitarity condition arising 
from intermediate states with more than two pseudo-Goldstone particles only start at 0{p^) — cf. [66l - [68| .) These 
amplitudes are constructed along the same lines. The details of the construction, including a full isospin breaking 
arising from 111^^0 ^ m^± , will be published elsewhere [36j (cf. also (STjV In this work we concentrate on the qualitative 
description of the result. 

The dispersively constructed scattering amplitudes of the pseudo-Goldstone bosons (pGB) take the following general 
form 

Ais, t, u) = J\f {V{s, t, u) -f U{s, t, u)) + 0(/). (65) 

Here M is an overall normalization and V{s, t, u) is a third order polynomial with the same symmetry properties with 
respect to s, t and u as the complete amplitude A{s,t,u). The coefficients of these polynomials for the independent 
amplitudes related by two-particle unitarity in all the crossed channels are identical for all the amplitudes and are 
the only free parameters entering the game. The non-analytic unitarity part U{s,t,u), which takes into account the 
contribution of the two-particle intermediate states in all the crossed channels, is then a known function of these 
parameters. In the low-energy region, intermediate states containing more than two pGB states contribute only to 
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the 0{p^) remainder, while intermediate states involving other hadronic states contribute to the coefficients of the 
subtraction polynomial. 

In the case of the amplitudes concerning one rj and three pion states, there are several two-pGB intermediate 
states to consider: tttt, KK, nrj. Since we shall only be concerned by the decay region, only the nearest singularity, 
coming from the cut produced by the tttt intermediate state, will be close enough to affect sizably the amplitude. The 
contribution from the remaining states {KK, nrj) can be expanded in a polynomial, which is included in V{s, t, u) (see 
also the discussion at the beginning of the next section). Of course, such an approximation would not be appropriateF^ 
to describe the rjir — tttt amplitude in the scattering region. In conclusion, for our purposes the only relevant related 
amplitude is therefore the tttt scattering one. 

For the charged 77 — >■ Stt decay channel the polynomial Vxis, t, u) can be expressed in terms of six free parameters 
corresponding to the t — u symmetric expansion at the center of the Dalitz plot 

Vx{s,t,u) ^ Axml+Bx{s-s')+CAs-s'f+Dx [it - s'f + {u - s'''f]+Ex{s~s''f +Fx [{t - s'f + {u ~ s'f] , (66) 

which is closely related to the traditional PDG parametrization of the Dalitz plot distribution. We take the overall 
normalization as 

so we have simply (cf. ([T91) ') 

Mx{s,t,u) ^Vx{s,t,u) + Ux{s,t,u). (68) 

Let us make one remark concerning the Dalitz plot parametrization. Between the polynomial (|66|) and the linear 
parametrization (|24l) there is a simple connection. However, the dependence of a, 6, d, /, g on parameters A^, B^, Cx 
and Dx is complicated by the presence of these four parameters also in the unitarity part Ux{s, t, u) (see below). The 
direct correspondence can be, however, established for the dependence of /, g on Ex and Fx with very simple form 

f~g^Ex + 2Fx, (69) 

which we will need in Sec. IVI Al (the exact connection will not be needed). 

For the related tt+tt^ — >■ 7r°7r° scattering amplitude (which is the only independent one in the isospin conservation 
case) we choose the following parametrization of the polynomial part in terms of the subthreshold parameters [66l \q1\ 

V..{s, t, u) = ^a^Ml ^A-A +j^{s~ 2M^f + |l [(* ^ ^AfJ)' + {u - 2Mlf 



{t~2M^^)' + {u~2^Q' (70) 



and the overall normalization A/Vu- — Ftt'^- The unitarity part of the 77 — > Stt decay amplitude Ux is then a function 
of a subset of the above polynomial parameters, namely 

Ux ^Ux {Ax, Bx,Cx,Dx;a^, 13^, Xi,X2). (71) 
The general form oiU{s,t,u) for the process AB CD reads 

UAB^Gii{s,t,u) = VK^(s) + W^{t) + W^{u) + {t- u)Wl{s) + {s- u)W^{i) + {t - s)WIj{u), (72) 

where the discontinuities of the functions Wg'}pjj{s) are given in terms of the right-hand cut discontinuities of the S 
and the P partial waves Si, Tg and Ui {£ — 0, 1) of the processes in the s-, the t- and the it-channels, respectively, 
aS 

disc W^{s) = 16ttI disc So{s)+ Aab Acd disc /'^""^^ ) , (73) 
V Vb(s)Acd(s)/ 

disc W^{s)= 48tts disc ^/^^^''^^^ , (74) 



However, the presented construction can be extended also to include the unitarity cuts from the other two-pGB intermediate states which 
are relevant in the scattering region [6^ . This then, however, brings into a game more free parameters (describing such intermediate 
processes). 

Note that in our case we need to continue these discontinuities analytically and they become complex (cf . [27l |3^ . |69| . [70| ) . 
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where and Aij(s) were defined in (|16|) . Similar relations for W^'^(s) can be obtained by an appropriate permu- 
tations of A, . . . ,D. The right hand cut discontinuities are iteratively constructed from the generalized two-particle 
partial- wave unitarity relations as described in [sB-Ell. The functions Wg'];^^{s) are then reconstructed by means of 
appropriately subtracted dispersion relation. Note that such a subtraction prescription is an indivisible part of the 
definition of the polynomial part of the amplitude. The first iteration reconstruct the amplitude at 0{p'^) while the 
second one yields the 0{p^) results. 

For the rj — > Tr+Tr^Tr*^ decay the above general form simplifies since there are only two independent masses in the 
problem and the amplitude is t — it symmetric. We get 

U^{s,t,u) = VF^(s) + W^{t) W^{u) (s - u)W^{t) + (s - t)W^{u), (75) 

where the subscripts S*, T refer to the ryvr" — > tt+tt^ and the 7/7r+ — > tt+tt" channels, respectively. The relevant 
discontinuities can be rewritten schematically as 

5 4 



discW^°(s) = ^J-.(.) 



fe=-l 



and similarly for (with coefficients a^*''), while 



i=l \k=-l ^ ' k=-l ^ ' k=G 



(76) 



(77) 



Here o^*\ a^*\ . . . , d!"^^ are known polynomials of the parameters {A^, B^, Cx, D^] ajr, /^tt, Ai, A2} and the masses m^^, 
m^; A(s) = A^^(s), a{s) = Xtt-k{s)/ s and Fi(s) represents a set of elementary functions listed in Appendix A. 

The corresponding functions Wgj,jj{s) are now expressed in terms of the dispersion integrals (the Hilbert trans- 



forms) Qi{s) of these functions, i.e. 



^' f°° dx J"^(a 



(78) 



with a suitable number ki of subtractions, < fci < 1. The 5- wave contributions in the s- and the t-channels are 
given by 

W|{s) = j2G^{s) J2 (79) 



5 



4 



i=l k=-l 

The P-wave contribution in the t-channel is more complicated. 



(80) 



i=l 



fe=0 



s) cfs^ + 



Y§l'\s)c^h^ (81) 



where (in the following formulae to± = to,, ± m-rt) 



,(A) 



r(A) 



— TO^j_ s — m_ 

g^{s)-g,{AMi) 

s - 4M2 

g\^\s)~g\^\Q), 

1 1 / 2 ^i(fli_^i(!^ _ 2 g»(g)-g»(™-) 

s — TO::j_ s — m_ 



g^(s) 



9 9 
TOq_ — TO_ 



to: 



g^is)~g,{ml) 4 g^{s) - g,{mi) 



rn_ 



s — TOl 



s — m_ 



(82) 

(83) 
(84) 
(85) 

(86) 
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The dependence of these functions on Gi{s) ensures the correct discontinuity of the function W^{s) and in addition 
is dictated by the requirement that the appropriate behavior in the chiral hmit [s^, |63] is reproduced. 

The expHcit form of the functions Gi{s) as well as the properties of the Hilbert transform are discussed in Appendices 
A and B. Here we only illustrate the above general procedure by means of the explicit result of the first iteration 
corresponding to the O(p^) part of the amplitude, and briefly discuss the 0{p^) result. 

A. rj 3n at one-loop order 

At the one-loop order our dispersive representation (j68|) of the amplitude A4x{s,t,u) simplifies substantially. The 
polynomial Vx{s, t, u) is only of the second order, 

t, u) = A.rn^ + - s') + C^{s~ s'f + [(t - s'f + {u ~ s'f] , (87) 

and of all the functions Ti{s) and their Hilbert transforms Qi{s) that were introduced in the previous section only the 
case i = 1 occurs in the unitarity part (j75p . Besides we only need the first term from (|8ip . 
The single function appearing at 0{p^) is thus 

TT 74JV/2 X X — S TT 

The form of Fi was chosen in order to ensure the relation Qi = (IBtt) J(s) [g^l (known also as Chew-Mandelstam 
function (tH]). 

The form of the unitarity part (j75p at the 0{p^) order is extremely simple in this formalism. For the polynomials 



introduced in (|80)) and (|8T|) in the case of the charged decay -H- tt+tt 7r° we find 

IGTrF^^^' - \A.,{7a^ - l&P,)mlml - ^B,(a^ - I3,)ml[iml + ml), (89) 
D y 

leTrF^a^i' = 2A.,p^ml + ^B.,{Aa^ml ~ P^{7ml + m^)), (90) 

167rF>W = ii?,/3^ (91) 

for the polynomials of S-wave in s channel (the polynomials that are not displayed here are identically zero). And 
then similarly for S-wave in i-channel 

167rF24'' = ^A,{a^ + 2/3^)m^m2 + ^B,{a^ + 2/?^)m2 (Sm^ + mf,), (92) 

ISTrF^a^'' = -]^A,(3^ml - ^B.,{2a^ml + p^{7ml + m^^)), (93) 

IQirFlali^ = is,/?, (94) 
and finally the polynomials of the P-wave contributions that are not zero are given by 

le^F^^foW = -^-B^ml, (95) 

le^F^^foW ^ -1-B,J3^. (96) 



B. — >^ Stt at two-loop order 

The whole amplitude at two loops, or equivalently at Oijfi) order, is of course more complicated. We employ here 
the full form of the polynomial (|66l) . The non-trivial part follows from the same general form (175^ with the functions 
Wgj.jj{s) from (|M 1) - ([5T|) . but contrary to the one-loop situation, where we have only one function Qi{s), we have 
to deal with five basic functions Qi, together with five derived types ([5^ - (l551) . Let us explicitly write down the first 
coefficient (which stands in front of Qi in the s-channel of S partial wave and thus together with (j89p represents the 

full a^^^ at 0(/)): 



a{s) loj 



<s) ~ 1 
ct(s) + 1 
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(16^)3a4^) = ^^A^^ml [(1152^2 - 2619a,/3, + SISO/?^) - 45/3,(a, + 2/3,)m2 + 11527r2 (17Ai + I8A2) ml] 



27 

81Fi" 



+ -j^B.m^ [/3,(61a. - 514;3^)m^ + (-459a^ + 1170a^/3^ - 1996/3^) m^m^ 



3 (447a2 - 987a^/3^ + 2150/3^) - 460871^ (2Ai + 3A2) ml {iml + m^)] 

2LFf' 
256 
27^?' 



+ [7a- (3"4 + T^IY - 2/3- (81"4 + SOrn^m^ + 17™;;)] 

+ ^^^""D^ml [7a^ {3ml + m^f - (3^ [lllml + 42mlml + 43m^)] 



(97) 



(i) -ii) 

From this example one can infer the general structure of all other parameters aj. , . . . ,df. . The full form can be 
obtained from the authors upon request. 



V. CONNECTION WITH CHPT: ORDER-BY-ORDER CORRESPONDENCE 



Let us briefly comment on the connection of the dispersive construction with the standard ChPT expansion. In 
analogy to the dispersive one, the 0{p^) ChPT amplitude can also be split into a polynomial part and a non-analytic 
unitarity part. The former corresponds to the tree- level counterterm contributions as well as to the chiral logs and 
sunset graphs, while the latter takes explicitly into account the nontrivial contributions of the loops. Though this 
splitting is not unambiguous and depends on the particular definition of the nontrivial part of the loop graphs, the 
unitarity part has to reproduce the correct discontinuities of the amplitude as required by (generalized) unitarity and 
corresponding to the two-particle intermediate states. Along with the pure pion loop contributions also the higher 
intermediate states are taken into account, namely, the graphs with kaons and rj inside the loops. However, below 
the 7777 threshold the contributions of discontinuities corresponding to the Trrj, KK and 7777 intermediate states are 
analytic and can therefore be expanded in powers of s, t, u. Sufficiently far below these thresholds one can show 
that their effects can be approximated by means of only the terms up to the third order (cf. [t^ and the numerical 
estimate of such error made in Sec. lVI A| . As a result we should obtain in this region an approximate ChPT amplitude 
with the same structure as our dispersively constructed amplitude (recall that both of them include the higher non- 
Goldstone intermediate states contributions only effectively through the low-energy and the subtraction constants, 
respectively). The only difference is that the polynomial part of the 0{p^) ChPT amplitude is generally complex 
due to the contribution of the sunset diagranrflwith three intermediate pions which develop nonzero imaginary part. 
However, it has been found to be tiny in |2C)I. l72j and therefore can be neglected. We reverify this observation in 
Sec. lvra 

These common features of both amplitudes suggest that the 0{p^) ChPT amplitude .4chPT(s, t, u), which we write 
in the form 

AchPT{s,t,u) = Af^(^M''clpT{s,t,u) + M''clpT{s,t,u) + m'^c^^ (98) 

can be reproduced as a special case of the dispersively constructed one. This can be quantified as follows in terms 
of what we call order-by- order fit. The ChPT amplitude in our dispersive parametrization is then represented by 
expressing particular chiral orders of our subtraction constants A^, . . . ,Fx and a^, . . . , A2 in terms of the LECs of 
ChPT, quark masses and chiral logarithms. Such expressions are then useful when one wants to organize the chiral 
result and to identify the renormalization-scale invariant combinations of LECs on which the amplitude depends. For 
the aims of the current work, it is however sufficient to perform this matching numerically and obtain the numerical 
values of our subtraction constants using the procedure described in the following lines (note that the same procedure 
would remain valid also if we wanted to obtain the analytic expressions, but instead of fitting the numerical results 
we would just compare expressions coming from ChPT with the ones of the analytic dispersive construction). 



Note that this diagram (77 — > Stt — > Stt) does not contribute to the unitarity cut of the riir — nn amplitude but instead its contribution in 
the decay region is analytic and can be expanded into polynomial. This polynomial can be complex since tti^ is unstable (tti^ > 3m.,r). 
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Let us formally split the parameters A^, ■ ■ ■ ,Fx of our amplitude into their 0(p^), 0{p'^) and 0{p^) parts, i.e. 

(99) 
(100) 
(101) 
(102) 
(103) 
(104) 

This induces a following splitting of the polynomial part of the amplitude 

nAs, t, u) = (s, t, u) + V^^^ (s, t, u) + T^^s) (s, t, u), (105) 



A 


- ^(2) , 
^x ' 




- aa(^) 




— ^x ^ 




- AS^*^) 


Cx 








Dx 








Ex 








Fx 









where 



V'^\s,t,u) = A'^i^ml + B'i\s - s^), (106) 

7'(^)(.,t,^) = A44)™2 + AB(4)(^ _ ^c) ^ ^(4)(^ _ ^,c)2 ^ ^(4) [(^ _ ^,c)2 ^ _ ^c)2] ^ (^(37) 

+ 4«)(s - s^f + F^'^ [{t - s^f + {u- s^f] . 



(108) 



Note that the unitarity part Ux{s,t,u) splits by construction naturally into the genuine one-loop O(p^) and the 
remaining 0{p^) parts that correspond to the first and the second iteration of the generalized unitarity relations, 
respectively, (see [30| for more details), 

= Ui^^ iA,,B,;a^, (3^) +ll<f^ (A,, B„ C„ D,; a,, Ai, A2) . (109) 

The unitarity part ui'^"' has been given in Sec. IIV A[ where we have written out the explicit dependence on the 
polynomial parameters of the 77 — >■ Stt and tttt tttt amplitudes. The 0{p^) part consists further of the genuine 
two-loop part and the one-loop part 

U^'^\A,,B,,C,,D,; a^,(3^, Ai, A2) = U^%^^^{A,, B,; a^, /?,) +u[%^^^{A,, B,, a„D,;a^, /3,, Ai, A2). (110) 
The 0(^2) ChPT amplitude M'clpr^{s,t,u) is now exactly reproduced by Vx'^\s, t, u) with 

4^) = I^i^, 5(2) = 1. (Ill) 

The imaginary part of the 0{p'^) ChPT amplitude below the irrj threshold is fixed by unitarity and therefore there 
holds exactly 

lmM^^lM^,t,u)=lmUi''HAi'\Bi'^;a<^\l3i'^), (112) 

where 

= /3(2) = 1 (113) 

are the leading order ChPT values of the tttt — tttt subthreshold parameters. Hence, up to a polynomial of the 

second order in s, t and u, the amplitudes A^Qi|prp(s, t, u) and (^i^^ , B^'' ; ai^^ , p'i'^) coincide (here we have tacitly 

assumed that the higher two-particle intermediate states contributions to A^|^]|prp(s, w) has been expanded to the 
second order in s, t and u as described above) and we can therefore write 

M'-^l^^is, t, u) = T'^^) (s, t, u) + Ui^^ (4'^ , Bi'^ ; a^^) , ) (114) 
for appropriate AA^x^ , AB^^^ , ci^) and ui^' . These parameters are found numerically by fitting the difference 

A(4) (., t, u) = M^^lM^, t, u) - Ui^^ (42) , B(2) . 42) ^ ^(2) ) (^^5) 
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to the second order polynomial Vx {s,t,u). When these 0{p'^) parameters are fixed, we proceed similarly to the 
0{p^) order. We compute the 0{p^) corrections to the unitarity part, 

where in addition to the parameters known from the previous steps there appear the NLO corrections of the subthresh- 
old parameters of tttt — )■ tttt scattering that are needed as inputs to this procedure. The discontinuities originating 
from the tttt intermediate states in s-, t- and u-channels of Ai^lp j,{s,t,u) and of this Vx^\s,t,u) coincide (modulo 
a power expansion of the higher-intermediate-state contributions to the third power as discussed above) . Finally, we 
fit the difference 

A(6) (s, t, u) = M^^lp^is, t, u) - V(6)(s, t, u) (117) 

to the third order polynomial Vx^\s,t,u) and set the remaining 0{p^) parameters. In this way, all the parameters 
of the polynomial part of the amplitude are numerically determined and the 0{p^) ChPT amplitude y^chPT(s, t, u) is 
represented now as y^chPT(sj t, u) — > .4|^Jd''(s, t, u), where 

At'^{s,t,u)^Vi'Hs,t,u)+Vi^H^,t,u)+Vi'Hs,t,u)+Ui''\A(^\Bi^^;a^ (118) 

By construction, the chiral orders of the various contributions to AQ^r^{s,t,u) were strictly respected — for 

(2) (2) 

instance the genuine two- loop unitarity corrections depend only on the leading order parameters Ax , Bx and 

(2) (2) 

air , /^TT • However, the known general form of the dispersively constructed amplitude A{s,t,u) can be further used 
in order to go beyond the strict chiral expansion and partially resum also the higher chiral-order contributions. This 
representation that we call resummed fit can be achieved by means of inserting the full parameters Ax, ■ ■ ■ ,Fx obtained 
by the above order- by- order fit and the full 0{p'^) tttt subthreshold parameters (or even the experimental values of 
the TTTT subthreshold parameters from [tI]) into the unitarity part of the amplitude, i.e. to define 

Af^lPis,t,u) ^Vi^\s,t,u) +ri''Hs,t,u) +ri''\s,t,u) +Ui^\Ax,Bx;a^,M +Ui''\A 

(119) 

The difference ^J?cs^~-^oiT order 0{p^) and contains effectively contributions of the one and the two- loop graphs 
with higher-order counterterms. It might be therefore treated as a rough estimate of the convergence of the chiral 
expansion. 

Let us note that we could also use another parametrization of the relevant tttt scattering amplitude based on the 
scattering lengths and effective ranges instead of the subthreshold parameters (see ^36, 37] for details) and repeat the 
above construction along the same lines. In such a case the amplitude Ajj^J{s,t,u) has to be numerically the same 
as before, namely, the parameters Ax, . . . ,Fx should be the same. However, the amplitude will now depend on the 
scattering lengths and the effective ranges of the tttt scattering taken up to the order 0{p*). Provided we then use 
the experimental values of these parameters in the resummed amplitude Af^^^{s,t,u), we can interpret the result as 
a partial resummation of the two-particle rescattering in the final state. The numerical effect of such a resummation 
might be even larger than within the previous parametrization, because the scattering lengths are known to have 
much worse convergent chiral expansion than the subthreshold parameters. 

VI. ANALYSIS OF THE CHARGED DECAY: ^ tt+tt tt" 

We have prepared everything to employ the dispersive representation for our analysis of the process rj — > Sir. It 
proceeds as follows. We start with the NNLO result of ChPT We determine the values of our parameters 

that reproduce the ChPT result, thereby checking also that the correspondence between these two frameworks holds 
using the order-by- order fit as outlined in the previous section. Our further analysis is motivated by the conclusion 
of Sec. nil Al that the observed mismatch between the ChPT 0{p^) predictions of the Dalitz parameters and their 
experimental determination by KLOE might be caused by the incorrect determination of the 0{p^) LECs Ci of ChPT. 
We therefore study the dispersive representation of ChPT with the values of the C^s undetermined and try to find 
the values of their combinations that reproduces the experimental data. Finally, after that we change completely 
the strategy and fit directly our dispersive representation to the experimental data. Such a fit gives us the 77 — )■ Stt 
amplitude up to the normalization that is determined from the matching with ChPT in the region where we can 
believe the ChPT result. In all the cases we are interested in the distribution we obtain and then by comparing the 
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FIG. 1. Domain of applicability of the dispersive method together with the physical region (ellipse). Full magenta line represents 
the nrj threshold in the M-channel (the only tttj threshold visible for this range of s and t). The blue dot-dashed line represents 
the axis of the t — u symmetry of the amplitude: e.g. the Adler zero in s, t = (represented by a circle) has its counterpart for 
s = and t = + 3m^. Dotted lines denote the 4m^ thresholds in the s- and the t-channels. Dashed lines define different 
regions in the matching procedure (see main text). 

decay widths computed from these distributions (by integration of the square of the amplitude over the physical phase 
space) with the experimentally measured one, we obtain the value of R. 

In principle, this could be done for both the charged and the neutral 77 — )■ Stt decay. However, as was discussed in 
Sec. IIIB[ no current experiment determined more than just one Dalitz parameter describing the neutral decay, thus 
we concentrate mainly on the charged one. Even in the charged sector the experimental situation is poor — only 
KLOE '22'] published just 4-1-1 Dalitz parameters (the last one claimed to be compatible with zero) describing the 
amplitude. From these values of the Dalitz parameters we have constructed a distribution in the physical region (in 
the similar way as done in (t^ ]) and all our experimental fits arc fits to such KLOE- like distributions, in our analysis 
we therefore depend fully on these KLOE measurements. 

A. Order-by-order correspondence: obtaining numerical ChPT distribution 

As was discussed in the previous section we can obtain the approximate chiral 0{p^) amplitude as a special case 
of our dispersive parametrization with some particular values of our parameters. The correspondence between such 
amplitude and the result of ChPT has to be almost identical neglecting only small effects descending from expansion 
of the two-kaon and the ttt] contributions and a tiny imaginary part produced by sunset-like diagrams. In principle, 
working in the (s,t) plane they should agree in the regioiQ for small s,t under the nrj thresholds in all the crossed 
channels. Although much bigger deviation should be visible only after KK threshold (the contribution of ttt] is very 
small) we stick on this as a strict limit of our method. Influence of systematic uncertainties is studied using different 
regions in our matching procedure (see below). The physical and the matching regions together with the nrj threshold 
are depicted in Fig. [1] 

We match the amplitude along the lines of the previous section really order by order. The correspondence for LO 
and the imaginary part of NLO can be verified analytically, having Ax and Bi from (|111[) and a\r — f^T, — I- 
After that we have proceeded with the matching numerically. From the NLO real part of the amplitudes we have 
fitted the parameters A^^^ , si"*^ , C^^^ , ui^^ (in the notation of Sec. El A^^^ = A^^^ + AAi^^ , etc.). After that we have 
verified the matching of the imaginary NNLO amplitudes and finally from the real part of the NNLO amplitudes fitted 
the parameters A^x\ ■ ■ ■ ,Fx^'^ (again A^x'^ — Ai'*' + AA'^'' ; this superscript is used for the NNLO values of ^2;, . . . just 
to distinguish these values from the ones of the overall fit from Sec. IVI Cp . 



Note that even though we are talking about the expansion for small Mandelstam variables (e.g. s and t), it does not simply mean that 
the smaller these variables are, the better agreement between these theoretical frameworks we obtain. The amplitude depends on three 
kinematic variables s, t, u which are connected by relation mOI I. So keeping two of them small, the third one is shifted up by the m,, mass. 
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set 


set 1 


set 2 


set 3 


set 4 


4(4) 
■^x 


0.464 ± 0.017 


0.457 ± 0.040 


0.459 ±0.011 


0.452 ±0.016 


0.462 ± 0.010 


d(4) 


1.95 ±0.15 


1.90 ± 0.21 


1.91 ±0.05 


1.88 ±0.09 


1.93 ± 0.06 




-0.42 ±4.4 


-0.68 ±0.91 


-0.62 ±0.17 


-0.76 ±0.41 


-0.51 ±0.25 




1.07 ±3.8 


1.04 ±0.03 


1.04 ±0.01 


1.04 ±0.02 


1.04 ± 0.01 


The values of the 0{p^) dispersive parameters 
amphtude. 


(in appropriate powers of [Gel/]) 


corresponding to 




set 


set 1 


set 2 


set 3 


set 4 




0.577 ±0.013 


0.581 ± 0.003 


0.581 ± 0.023 


0.577 ±0.002 


0.583 ±0.011 




2.42 ±0.22 


2.460 ±0.012 


2.47 ±0.19 


2.44 ±0.01 


2.49 ±0.10 


'^X 


0.24 ±3.4 


0.30 ±0.11 


0.38 ± 1.7 


0.20 ±0.09 


0.55 ± 0.90 




1.55 ±2.4 


1.57 ±0.02 


1.58 ±0.03 


1.58 ±0.02 


1.58 ±0.02 




5± 149 


5.4 ±0.7 


5.6 ±4.7 


5.1 ±0.6 


6.1 ± 2.6 




-4 ±84 


-3.6 ±0.1 


-3.7 ±0.2 


-3.7 ±0.1 


-3.7 ±0.1 



TABLE VII. The values of the O(p^) dispersive parameters (in appropriate powers of [Gel/]) corresponding to matching to 
ChPT NNLO amplitude. 



Concerning the tttt part we follow closely the determination of its subthreshold parameters as established in 67 j 
For the particular values we have used (|113p for the leading order and set the NLO values to be 



v(4) 



1.044. 



L083, 

A^^^ = -1.43 X 10-^ , A^^' = 8.5 x 10"^ 



(120) 



The fits were performed in the following regional (all numbers in GeV^; cf. Figure [T]): 

• set 0: the physical region; 

• set 1: the square region around the Adler zero {s,t) G (—0.05,0.05); 

• set 2: the triangle region between the lines s — Arri^, t — Arn^ and the Trry threshold; 

• set 3: s € (—0.05,0.05), t between the Trrj threshold and the t ^ u axis; 

• set 4: s e (—0.1,0.1), t between the ttt] threshold and the t ^ u axis. 

Distances between the points in grids are constant in both the s and the t directions, and the approximate total 
number of them is the following: set 1: 300, set 2: 900, set 3: 1600 and set 4: 4400. Further, for the physical region 
(set 0) we have chosen the same points that were used in [20i] . i.e. 174 points, which is a very similar number to the 
KLOE's number of bins (154, cf. also discussion in [T^I)- The different regions with the different numbers of points 
were set in order to have systematic and statical errors under control. 

For the fits we have used MINUIT package with the weights of the individual points set to AO{p^)/2. Results for 
the NLO parameters are summarized in Table I VII and the ones for the NNLO parameters then in Table IVIIl The 
error bars quoted for the individual parameters are results of MINUIT. 

At the moment we have in hands the dispersively constructed amplitude (i.e. the analytic formula) which is nu- 
merically equivalent (or very close) to NNLO ChPT amplitude. We can verify the equivalence also by computing the 
decay width. 

Our dispersive representation was constructed in accordance with chiral perturbation theory and we have chosen 
similar normalization as used in ^0^ (with an extra factor 1/F^). We can thus compare directly a neat amplitude 
Aixis, t, u) with the isospin-breaking factor pulled out as defined in ([T^. The result of the integration of the amplitude 
square over the physical phase space is (cf. (6.7) in (20j): 



r(r; 



tt") ^ sin^ e X 2.68 MeV, [ChPT] 



(121) 



Since the amplitude is t ^ u symmetric, one can fit it only in the region below the t — u line. 
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where we have introduced 

sine-^. (122) 

Comparing this result with the experimental measurement for the decay rate ^ij we arrive at the value which exactly 
reproduces the one of [l^l (mind the typo in [20| ) 

R ^ 41.3. [ChPT] (123) 

We can use result [20| also for a numerical estimate of the error induced by a few approximations in our parame- 
trization we have made with respect to the ChPT computation. As was discussed in Sec. |Vl we have neglected the 
imaginary parts of our parameters (which are connected with the contribution of the sunset diagram) . In the physical 
region we have performed fits, in which we have allowed the 0{p^) parameters to be complex. We have found that 

the NNLO ChPT result is very well approximated by adding a constant imaginary term Im^i^'' — 0.080 ± 0.064. By 
neglecting this term in the computation of R we introduce an error of 0.1%. Similarly, we have neglected higher than 
third order polynomial terms in the expansion of KK and irr] contributions (in the decay region). We can estimate 
the corresponding error by addition of some higher-order terms into the polynomial. The symmetries dictate that the 
fourth order polynomial would contain terms Gx{s — s'^)^ + Hx{s — s'^)^{t — s'^){u — s'^) + Ix {it — s'^Y + [u — s'^)*). 
From the dimensional considerations, the contribution of KK intermediate states into these parameters should be 
~ w4 ,} p )2 ^ 10^ GeV~^ (and similarly for ttt]), whereas even if all of them were ~ lO'' GeV"^ the shift in the 

determined R would be 0.1%. Both of the errors are therefore negligible with respect to the other sources or error 
discussed in the following analyses. 



B. Correction to order-by-order fit: Correcting the ds in ChPT 



In the previous subsection we have constructed the dispersive amplitude reproducing ChPT in the region where our 
method is applicable. It is no surprise that if we fitted this dispersive representation to the Dalitz parametrization 
p3| as was done in ^2d\, we would obtain the same values of the Dalitz parameters as [131. In Sec. lIII Al we have found 
an indication that the discrepancy between so obtained values and the values measured by KLOE can be (at least 
partially) caused by the incorrect values used for the 0{p^) LECs Ci of ChPT. The contribution of the C^s to the 
0{p^) amplitude is polynomial and real and so changing them means changing the 0{p^) part of our polynomial (jlOSp 
— shifting the parameters appearing in it. By studying the chiral amplitude obtained from our previous analysis with 
an unknown 0{p^) polynomial added, 

A'x'^'is, t, u) = Axis, t, u) + AAxis, t, u), (124) 
with AAx{s,t,u) = J\fjjAVxis,t,u) and 

AVxis, t, u) = AAxM^ + ABx{s - s") + ACx{s - s'f + ADx [{t - s'^f + (u - s')^] 

+ AExis - s'^f + AFx [{t - s'^f + (u - s")^] , (125) 

we can thus study the impact of the corrected C^s on the chiral 0{p^) amplitude. 

Provided the dominant part of the discrepancy between the NNLO chiral result and the measured amplitude is 
hidden just in the incorrect determination of the CiS, the chiral 0{p^) amplitude with the correct set of the C^s, 
and thereby also the corrected amplitude Ax°^{s,t,u), should reproduce the physical data. Therefore by fitting the 
KLOE-like distribution, we should obtain the values of the dispersive parameters corresponding to the correct values 
of the CiS. By comparison of these values with the analytic expressions of these parameters in terms of the C^s, one 
could obtain approximate constraints that the correct values of the C^s should fulfill^. 

However, in the case the change of the C^s is insufHcient in order to obtain the physically measured amplitude and 
there still remains a big difference between the physical amplitude and the one obtained from NNLO ChPT, the fit 



At the current level these constraints could be formulated in terms of reproducing the measured Dalitz plot parameters. For every such 
parameter by using relations of Section III. A and the observed difference between its experimental value and the value coming from [20| 
with all d = 0, one obtains one constraint on the C'iS. Note that provided the information on R was supplied from another source with 
enough accuracy, one could obtain one additional constraint on the CiS. Unfortunately, such constraints are very complicated and would 
need to be analyzed together with additional constraints coming from other processes (similarly as was done in [65|) in order to provide 
any useful information on the values of CiS. 
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# 174 


# 2500 




-0.05 ±0.3 


-0.029 ± 0.003 


AS, 


-0.5 ±1 


-0.46 ±0.01 


Aa 


-7 ±2 


-6.97 ±0.07 


AD, 


-0.7 ±0.8 


-0.64 ±0.02 




-37± 18 


-36 ±3 


AF, 


24 ±5 


24 ± 1 



TABLE VIII. Corrections to the chiral dispersive parameters in order to reproduce the KLOE data (in appropriate powers of 
GeV). 

of the amplitude A'^'^ {s,t,u) to the physical data would mean this time that this difference was parametrized (and 
approximated) by the polynomial l!ii.Ax{s^t,u). 

Unfortunately, just from the fit of the amplitude y^^°''(s, t, u) to the physical data there is no way how to distinguish 
between these two scenarios — either that the 0{p^) chiral amplitude with the correct values of the Ci describes well the 
physical amplitude or that the higher-order remainder can be on the physical region approximated by the polynomial 
(or if both of the situations are present in a combination, there is no way how to separate these two contributions). 
We should be therefore careful with the interpretation of the result of such a fit and take this fit just as a starting 
point for the deeper analysis of the chiral amplitudes. Note that in Sec. IIII i^ we have listed a few criteria that would 
indicate the realization of the first scenario. 

In any case the amplitude constructed that way should describe the data better than the ChPT parametrization of 
the previous subsection and the determination of R from this distribution will be closer to the real one. 

We have therefore fitted all A parameters from (I125P using the condition that the amplitude (|124l) has to fulfill 
the distribution of data based on KLOE. This distribution is limited only within the ellipse of the physical region (cL 
Figure [T|). We have performed the fit for the following two data sets. First we have created exactly the same points as 
done for set (or equivalently in [23|) and then much more (2500) in order to study how this affects the dependence 
on the statistics. Using the different sets of parameters for the chiral amplitude Ax{s, t, u) as summarized in Table IVll 
and IVIII has very little effect on the resulting A parameters, so we display only their average, cf. Table IVIIII For 
further applications we have also fitted the complete 0{p^) polynomial v!f\s,t,u) + AVx{s,t,u), whose results are 
presented in the first column of Table HXl Since we have no further information from KLOE concerning the efficiency 
of every selected bin, the appointed errors in both tables are not very reliable. 

As we have discussed above, it is not easy to interpret the result (moreover, if the fit relies just on the KLOE-like 
distribution we have made just from the 4+1 measured Dalitz parameters and do not have any notion of the systematic 
uncertainties here), but let us present some interesting observations. All the corrections of the parameters seems to 
reduce the original 0{p^) contributions to these parameters coming from ^2^, some of them even tend to change the 
sign of the total 0{p^) contribution with respect to the original one. Taking into account relation together with 
the condition ([Bn]) imply that in the case that all the difference between the physical amplitude and the original 0{p^) 
chiral one is hidden in the wrong determination of the C^s, 

AEx + 2AFx = 0. (126) 

From Table IVIIII we see this tendency. 

Let us return to our main interest, the determination of the ration R. The integrated decay rate computed from 
this distribution is 

r(?7 ^ 7r+7r-7r°) = sin^ e x 2.24(10) MeV. (127) 

Together with the experimental input for F this leads to 

i? = 37.7(9). [ChPT±disp.±KLOE] (128) 

The quoted error is based only on getting values for parameters in different regions as explained in the previous text. 
It does not take into account systematic errors coming from the experimental data, which we do not know, and the 
theoretical error from the ChPT part. The later one can be estimated from the convergence of results coming from 
the chiral expansion. 

The LO value coming from the current algebra is i? = 19.1 and the NLO result was R = 31.8. Taking these values 
into account, our prediction for this ratio from the NNLO ChPT and data of KLOE is 



R = 37.7 ± 2.8. [ChPT+disp.+KLOE] 



(129) 
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cor. set 


fit to KLOE 




0.575 ± 0.006 


0.575 ±0.001 


Bx 


1.99 ±0.04 


2.15 ±0.02 


Cx 


-6.8 ±0.3 


-5.8 ±0.2 




0.94 ±0.03 


0.87 ± 0.08 




-31 ±3 


-19 ±9 




20 ± 1 


21 ±5 



TABLE IX. Dispersive parameters as free parameters fitted to KLOE (ultimate column; in appropriate powers of GeV). For 
a comparison in tlie first column we have also summarized the corrected parameters of the previous subsection. 

C. Overall fit: using the measured distribution 

So far we were using the dispersive representation in a very close connection with ChPT and up to the fact that we 
were trying to correct it by the values of the C^s coming from KLOE, we just reproduced the amplitude coming from 
this theoretical framework. But we can also change our strategy totally; we can use the experimentally measured 
distribution and employ ChPT just for the normalization. In this analysis we therefore assume that no matter what 
the proper description of this process leading to the correct physical amplitude would be, it would fulfill the general 
principles of quantum field theory together with the observed hierarchy of various contribution^^ which were used 
for the construction of our analytic dispersive representation. It then means that such correct physical amplitude can 
be to a good extent described by our parametrization and the later can be used for its analytic continuation to some 
region where ChPT gives a reliable result for the amplitude, and can be used there for the matching. Such procedure 
will substantially reduce the influence on R of the error connected with the chiral expansion of the amplitude. 

Let us start with the easier part — fitting the KLOE-like distribution, which plays in this analysis a role of 
the experimental distribution. In contrast to our previous fits, where respecting the chiral orders of the dispersive 
parameters was natural (and important), in this case keeping the different chiral orders of the parameters makes no 
sense. The more natural approach is using our representation in the resummed form — the values of the parameters 
in the polynomial and in the unitarity part are the same. 

The fit of this general representation to KLOE-like distribution was performed for the same data set as in the 
previous subsection (174 data points). Note that in the fit, the overall normalization is set so that the amplitude is 
equal to one at the center of the Dalitz plot. In order to simplify the comparison between this and the previous fits, 
we have decided to multiply all these data by the numerical factor that produces the same number for parameter 
Ax as the corrected value obtained from the analysis in the previous subsection (average over all data sets), i.e. 
Ax = A^x^ + AA, = 0.575. 

The values of so normalized dispersive parameters coming from the overall fit of KLOE are presented in the ultimate 
column of Table IIXI For the comparison there are displayed also the corrected values of these parameters from the 
previous subsection. Note that these two sets of values obtained from fit to KLOE correspond to two amplitudes with 
different unitarity parts (the first one contains the parameters of ChPT respecting their chiral orders, whereas the 
later contains in the unitarity part exactly those values of the parameters appearing in the polynomial part) and that 
in both cases the unitarity part forms an important part of the amplitude. The agreement between these two sets 
that is seen in this table is therefore quite interesting (the small disagreement is seen only for parameter which 
is, however, given with the biggest error in both approaches). Let us stress once more that the particular overall 
normalization was taken only to simplify this comparison and we still have to remember that values in Table IIXI are 
multiplied by, for the moment, unknown constant. 

To set this overall normalization is, in our opinion, the main issue of the dispersive study for rj — Stt, so we try to 
be as cautious as possible. First of all, we will rely on a set of points rather than only on one point even though it 
would be sufficient for setting the normalization. Thus, we need to select the region of the points where we believe 
ChPT result. In order to achieve this task, we discuss the following articles which could be important for its selection: 

i) Adler zero condition 

ii) correspondence with order-by-order fit and convergence of chiral orders 



This hierarchy is expressed in the construction of the representation in terms of the very basic chiral counting of the partial waves of the 
amplitude - cf. relations (2.2) and (2.3) in [3711 . 
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iii) plateau argument 



Let us explain them in detail. The point i) is connected with the SU{2) x SU{2) theorem which was summarized 
in Sec. Ill CI According to this theorem the values of the charged amplitude at the points in a small vicinity of 
SU{2) X SU{2) Adler zero, are O(m^) for rriTr — >■ and therefore protected from being large. At 0(p^) the zero of the 
real part of the amplitude, namely, the point s = u — lArn^ belonging to the above set of protected points, possesses 
the following additional convenient property. The slope of the amplitude (which is generally not protected by that 



theorem) develops O(p^) corrections that are accidentally small. In the previous dispersive analysis [28|, this fact was 
the main motivation for matching of the dispersively constructed amplitudes exactly at this point. However, as we 
have learned from Sec. Ill CI there is no guarantee that at the points near the Adler zero the chiral corrections to any 
calculated order are small. Nevertheless, the points where the real part of the amplitude vanishes (which are often 
called "Adler zeros" in this context too) can serve as good reference points, or benchmarks, of the individual chiral 
orders. In such a way these points were used in the analysis of the NNLO ChPT calculations ^] with a result that 
the best convergence of their positions is observed on the lin^H t — u. On the other hand, according to the same 
analysis, the point s = u = 1.4m^ does not seem to be particularly stable with respect to the O(p^) corrections. The 
matching of the overall normalization at the "Adler zeroes" has however the advantage that at these points just the 
imaginary part of the amplitude is matched and therefore the uncertainty corresponding to the not well known 0{p^) 
LECs is eliminated (or suppressed when we match in their vicinity). 

The second article ii) can help us to reformulate the previous conclusion in different words. The order-by-order fit 
should be by construction very similar to the chiral expansion. The convergence of this expansion was crucial in [2^ 
for setting the matching point in s — u around the 0{p'^) Adler zero (s — 1.4mJ). However, following the detailed 
analysis of [2^ we have to conclude that this choice becomes to be very dangerous for matching at 0{p^) (for s ~ u 
the imaginary part at NNLO is even of the opposite sign when compared with NLO). Much better convergence when 
coming from NLO to NNLO is seen for t = u, which we use in the following. 

Finally iii) reflects the stability of the points within the given region or cut. Studying some physical observable 
(for example the decay width) as a function of the matching point one would expect a plateau behavior in the correct 
region. 

Using the previous arguments we fit the normalization of subthreshold parameters for t = u cut (below the physical 
threshold) matching only the imaginary part of ChPT amplitude where we interpolate between both variants, the 
standard and the resummed one (cf. also Fig. [2] below). Within the straightforward analysis one obtains 

T{r] 7r+7r-7r°) = sin^ e x 2.25(40) MeV, (130) 

which corresponds to the value 

i? = 37.8 ±3.3. [disp.+KLOE] (131) 

This number depends more strongly on the data than it was in the case of the order-by-order fit of the previous 
subsection. This is the reason why we have dropped "ChPT" in its description even though one should remember that 
it enters the determination of this number through the normalization as explained above. The sources of errors are 
thus of the two types: the uncertainties connected with the experiment — the uncertainties of the experimental data 
we fitted and the accuracy of their parametrization by our dispersive representation; and the uncertainties connected 
with the normalization procedure — the error that is induced by the analytic continuation of the parametrization 
from the region where we have fitted the data to the region where we have matched with ChPT and finally, the error 
of the the determination of the values of the amplitudes in this region from ChPT, the error of the chiral expansion. 

From the way our parametrization was constructed and from the normalization procedure described above, it 
should be obvious that we have concentrated mainly on the reduction of the errors of the second type. Moreover, 
since we used only a distribution coming from 4-1-1 Dalitz parameters measured by KLOE without any information 
of the systematic uncertainties in different regions of the Dalitz plot, it should be obvious that these uncertainties 
prevail and the error quoted in (jl3ip corresponds solely to them. Provided we have more precise measurement of 
the distribution of the physical amplitude, we could quantify the uncertainty coming from fitting these data to our 
parametrization and a deeper analysis of the ones coming from the normalization procedure would be required. 

To conclude the study on the charged decay rj iT^Tr~Tr^ let us summarize all the analyses performed in this and 
the previous Sec. I VI Al and I VI 51 in one plot (Fig. [2]) focusing on the t = u line. 



'Note that these points do not belong to the set 1391 1. i.e. they are not close to the Adler zero in the strict sense of the SU{2) X SU{2) 
theorem. 
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FIG. 2. The real (blue) and the imaginary (red) part of the amplitude for t = u line. The dotted lines represent order-by-order 
fit (Set 3), the dashed ones stand for resummed fit constructed from it (with the values of all parameters equal to their O(p^) 
values from Set 3), and finally the solid lines refiect an overall fit corresponding to KLOE values, with normalization set to 
interpolate between the dotted and the dashed lines for the imaginary part. The vertical lines indicate the physical region. 



•q — !> 7r^7r~7r'' 6 8 

77 ^ 37r" 5 7 

r; — >■ 37r 6 9 



TABLE X. Numbers of free parameters of the dispersive parametrization that are needed to be determined in the various 
analysis. The first two lines describe the individual studies of the charged and the neutral 77 — > 37r decays, while the ultimate 
corresponds to the combined fit of both of them. 



VII. NEUTRAL DECAY: r) 37r° 

Before we present our results for the ry — > Stt*^ amplitude, let us shortly discuss the number of free parameters of 
our dispersive parametrization for the 77 — > Stt decays that has to be determined from a fit. 

The polynomial part of the charged amplitude was given in For the 77 — > 37r° decay the situation is simpler 

since 

Vo{s, t, u) = Aaml + Co [{s ~ s^f + {t - s^f + {u~ s'f] + Ea [(s - s'f + {t- s^f + {u- s^f] . (132) 

These two decays are related by the 2-particle unitarity and so in the unitarity part of the 77 — >■ tt^tt^tt^ decay 
there appear 2 parameters from the neutral decay, whereas in the unitarity part of the 77 — >■ Btt'^ decay there appear 4 
parameters from the charged one. When one takes the full isospin breaking into account there is no further connection 
between these two decays and one needs to determine all these parameters appearing in the considered amplitude. 
However, in the case we work in the leading order of the isospin breaking, relation (j20p bounds these two amplitudes 
together and all the parameters of the neutral decay can be expressed in terms of the charged parameters. The number 
of the parameters needed to be determined in the various studies in these two cases are given in Table [Xj 

From the table and the present status of information we have from experiment on these amplitudes (as summarized 
in Sec. IIIBI) . it is obvious why we were focusing just on the charged decay modes of ij. In that case we have six 
unknown parameters in our dispersive formula (in the case 77z^± = 77^^o) which could be saturated by five known 
Dalitz parameters of this decay. On the other hand, although the neutral decay is theoretically much simpler (having 
less parameters and there is no P-wave contribution to the unitarity part), so far only one Dalitz parameter (a) was 
measured for 77 — Stt*^. The procedure elaborated in the previous section will not be thus very reliable in this case. 

Working in the m^± — m^o we can obtain the values of the neutral dispersive parameters from the values of the 
charged ones (using (|20l) ). from which we can compute the neutral Dalitz parameters. 

Doing so we get a = —0.042 for the order-by-order correspondence (ChPT with an additional 0{p^) polynomial) 
and a = —0.047 for the overall fit (fit to KLOE with the normalization from ChPT). This two predictions can be put 
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FIG. 3. The absolute value (squared) of the amplitude rj Stt" for t — u line. The notation is the same as in the previous 
graph: the dotted line represents order-by-order fit (set 3), the dashed line denotes a resummed fit (with 0(p®) set 3 values), 
and the solid line stands for an overall fit corresponding to KLOE value. With the vertical lines we have demarcated the 
physical region. 



together to: 

a = -0.044 ±0.004. (133) 

In these error bars the systematic uncertainties from KLOE are not included. This error also does not take into 
account the isospin corrections to the relation (1^0]) . 

The corresponding study of KLOE (2^ based on the charged decay mode leads to the value a = — 0.038(3)(12), 
which is also higher (in absolute value) than are the most precise values on this quantity (see Table HIl e.g. MAMI-C: 
a = —0.032(3)). Note that the similar study done recently in [34] obtained the same prediction at the leading order 

in TTTT rescattering taking using KLOE data: a^^gprp = — 0.042(2) (g). Employing higher orders to this prediction 
leads to aNREFT ~ —0.06, the value even more off the current most precise experimental determination. However, 
these two values were really based on the KLOE data. The independent determination of ^34] for the parameter a is 
in much better agreement with the present most precise experimental value (see e.g. Table UlI)) . 

Therefore, we have verified a deviation in the neutral Dalitz parameter a obtained from the measurement of the 
charged Dalitz parameters (using isospin relation (j20l) ') and the one obtained from the direct measurement of the 
neutral decay. This discrepancy can be attributed to the fact that we do not have at disposal the exact KLOE 
distribution but only the simplified one. The second reason is connected with the following disadvantage of the 
dispersive representation when connected to the isospin simplification m^± — m^o . This limit forces us to use only 
one value for the two masses of pions. In the case of 77 — > 7r+7r~7r'^ we use an average value of the charged and the 
neutral one, whereas in the neutral case there is used naturally the neutral mass itself. In ChPT we are free to make 
this change without changing the other parameters of the amplitude (at least in principle LECs do not depend on 
the light quark masses). This is, however, not true for the polynomial parameters used in the dispersive approach, 
which do depend on the masses. Indeed, repeating the calculation of a done for the "average" pion mass we obtain 
the value of a that is reduced by approximately 10% with respect to our result (jl33p . i.e. it goes in a good direction 
to the independent measurements. 

Similarly as in the previous section for the charged decay, we summarize our analyses for the neutral decay in one 
plot (Fig. [3]) again on the line t = u, now for absolute value squared of the amplitude. Note the only slight change of 
the slope (a) when passing from the order-by-order fit (or equivalently NNLO ChPT) to the resummed one (which 
is equivalent to a different resummation of the 0{p^) chiral orders in ChPT). A bigger change in a, connected even 
with a flip of the sign, is seen only after inclusion of the KLOE data. 

Before closing our study of the neutral decay mode let us return back to the resummed correspondence introduced in 
Sec. lVIC] The major problem addressed there, the overall normalization, can be totally ignored for a special quantity, 
the ratio of decay rates 

I (77 — )■ TT+TT 

where such normalization simply drops out. Using the numbers from fit to KLOE in Table HXl one obtains 

= 1.475 ±0.015, (135) 
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the value that is in an excellent agreement with the most precise measurement of [75|: rrj = 1.46(3)(9). We have verified 
that our result is stable against isospin corrections and is very close to the original ChPT value [20] r^^^'^ = 1.47. 
Note that the prediction of NREFT ^ is somewhat smaller = 1.40(1)(4). 



VIII. CONCLUSION 



The main purpose of this paper was a demonstration of the possibilities of using the dispersive approach in the 
case of the 77 — >■ Stt decays. In the advent of the new precise measurements which could be sensitive to higher-order 
isospin-breaking effects it is a reliable theoretical tool, based only on general assumptions, as relativistic invariance, 
unitarity, analyticity and crossing symmetry together with chiral counting, which ensures that our amplitude is valid 
up to and including 0{p^). It easily reproduces the ChPT amplitudes, providing them in a simple and compact 
analytic form. It is therefore useful for studies of the structure of the amplitudes and enables us to identify the basic 
independent combinations of the LECs appearing in them. In addition, it is more general and can be analytically 
continued outside the kinematic decay region. Further, it incorporates naturally the isospin corrections induced by 
the m^± — TOttO difference. 

However, in the present situation when experiments are limited to 4+1 (one parameter is zero or close to zero) 
Dalitz parameters for rj — 7r"'"7r~7r'^ and 1 parameter for rj — Stt'^ decay it seems reasonable to use only leading order of 
the isospin breaking. This leading order is hidden in the parameter R (see p9)) ) and Ai{s, t, u) is thus effectively taken 
in the isospin limit. Even such simplified analysis has an important theoretical outcome since by comparison of the 
experimentally measured decay rates of these processes with the predictions coming from theory, we can determin^H 
the parameter R, thereby obtaining information on the individual masses of to„ and m^. 

Such analysis was performed within NNLO ChPT, but the necessity of determination of the large amount of NNLO 
LECs of ChPT, together with the observed discrepancies between the values of the Dalitz parameters describing 
the energy dependence of the amplitude predicted by ChPT and those measured by experiments question also the 
accuracy of the result obtained for R. 

On the analysis of the charged data from KLOE we have presented two methods how to make use both the 
information we have from the NNLO ChPT and the one from experiment together in order to determine the corrected 
value of this parameter. The first one is motivated by the possibility that all the discrepancy between the NNLO 
ChPT predictions and the KLOE measurement comes from the incorrect determination of the 0{p^) LECs of ChPT 
(or more generally that all the discrepancy can be approximated well by a polynomial of the third order). Using it 
we have obtained the value R = 37.7 ± 2.9. The second method is based on the fit to the experimental data, in which 
the error of a slower convergence of the chiral amplitude is reduced by using ChPT just for the normalization of the 
parametrization in the region where the chiral expansion seems to converge well. Here the main source of uncertainties 
was due to the experimental error and we obtained R = 37.8 ± 3.3. Even though both of these results depend on 
KLOE, the nature of the dominant error in each of the methods is different and so we can combine them to obtain 
our final predictiorF"! 

i? = 37.7 ±2.2. (136) 

This value agrees very well with the lattice average prediction [121] and is compatible with the Dashen's limit at 
NNLO order. (However, it would correspond to large Dashen violation if one worked at NLO as studied in [6^ — cf. 
also Table 6 in [2(1]). Using relation (U) and the value of r = nis/rh = 27 A ± 0.4 provided by lattice simulation p^. 
we obtain 

g = 23.1 ±0.7. (137) 

Note that [77| and several non- lattice determinations (for the recent overview see [t^ [t^] ) point towards a smaller 
value for r, which will mean also a smaller value of Q, e.g. the value of r ~ 24 leads to Q ~ 21.5. For the sake of 
completeness, let us also mention the value Q — 20.7 ± 1.2 obtained from large deviation from Dashen's limit in [13] 



(cf. also 181|)- 



'Note that instead of R we could use in our normalization of the amplitudes the parameter Q as well. Such choice is favored in the 
analyses based on ChPT at NLO since at this order Q depends only on QCD meson masses 0] and is reasonably stable with respect to 
the Kaplan-Manohar transformation l76l of quark masses of ChPT. However, both these advantages of Q are lost when one includes chiral 
two-loop effects because the relation between Q and the meson masses gains noticeable r = iris/rh dependent chiral corrections at NNLO 
(cf. [63|). When matching the amplitudes with the results of NNLO ChPT .20j, it is more natural to employ the normalization containing 
R and assume that the Kaplan-Manohar ambiguity is fixed by the values of LECs used in that computation (for instance the value of 
Lg = stemming from large Nc considerations). 

'Let us emphasize once more that this number stands and falls with the assumption that the genuine physical amplitude is described well 
by the KLOE-like distribution constructed from the values of 4-1-1 Dalitz plot parameters presented by KLOE |22)I . It is therefore desirable 
to remeasure the charged r; — > Stt decay and repeat the performed analysis with the access to the real data. 
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In the case we want to fully complete the task set in the Introduction and determine the current masses of to„ 
and niii at some scale /i, we need as inputs the values of rh and rrig at the same scale obtained from somewhere 
else. Introducing the most recent averaged values of PDG [4l|, in MS scheme at the running scale /i — 2 GeV, 
nis = (100.6 ± 2.1) MeV and m ~ (3.8 ± 0.1) MeV into the definitions ^ and ^ as inputs, one obtains the current 
masses at this scale 

TO„(MS, n = 2 GcV) = (2.52 ± 0.13) MeV, 

(138) 

TOd(MS, ^ = 2 GeV) = (5.08 ± 0.13) MeV. 

Employing instead the averaged lattice results [l^l rUg = (94 ± 3) MeV and m — (3.43 ± 0.11) MeV leads to 

TO„(MS, (1^2 GeV) = (2.23 ± 0.14) MeV, 

TOd(MS, ^ = 2 GeV) = (4.63 ± 0.14) MeV. ^^^^^ 

Our results are compatible both with the averaged values of these individual masses quoted in PDG [ill and with the 
averaged lattice results of [l^]. Note that these two sets for each mass are correlated since the PDG values contains 
among others also the lattice results. 

We have also made a quick look onto neutral decay mode rj — > 37r'^. Assuming that the original KLOE measurement 
of charged Dalitz parameters a, 6, d, /, g is correct we have verified some deviation in neutral Dalitz parameter a 
obtained from our representation using them (and the isospin relation pO|)) and from its direct measurements. In the 
previous section we have discussed the possible explanations of this discrepancy. Note that the neutral decay rj — > 37r'^ 
hides more — in the physical region there occurs an intrinsic cusp, however, its position is naturally very close to the 
edge of the phase-space. Its shape is thus very suppressed and the techniques similar to the extraction of Kl — 37r° 
have to be employed. An unambiguous description of the amplitude including corrections induced by m^± — m^o 
difference can thus naturally help in this extraction. 

At this point let us mention also the second neutral Dalitz parameter /?, whose measurement is planned also in 
the forthcoming experiments (either WAS A or KLOE-II). Its theoretical determination in ChPT is not influenced by 
0{p^) LECs and is thus important for its consistency check. In Sec. IIII Al we also list the other combinations of the 
Dalitz parameters that are suitable for such a check and their evaluation can shed light on the discrepancies in the 
determinations of the Dalitz parameters — not only the discrepancy in a whose solution was proposed in [s^l but 
also the discrepancy in b that even in [3^ 1 remains. 

The last physical quantity we have discussed was r^, the ratio of the neutral and charge decay width, with the 
result 

= 1.475 ±0.015. (140) 

Our prediction, or more precisely the prediction of KLOE based on our dispersive parametrization, agrees not only 
with the most recent experiment but also with the NNLO ChPT calculation. Naively, seeing Figures [5] and [3] and 
the change in the absolute value for both amplitudes, one would expect some impact, but apparently these two shifts 
canceled out in the ratio. 
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Appendix A: Kinematic functions appearing in the dispersive representation 



Five basic functions Gi{s), i — 1,...5 entering general formulae ([5n|) - ((55)) for Wg rp{s) and W^{s) are given for 
complex s by the dispersion integrals with an appropriate number of subtractions ki 

y^[S) = / — , Al) 
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where !Fi is the set of generally complex functions describing the discontinuities of the amplitudes listed below. 
Therefore the functions Qi are analytic in the cut complex plane with cuts along (4toJ, oo) and with the discontinuities 
across these cuts 

disc a, (s) = ^ {g^{s + iO) ~g^{s~ lO)) = F^{x). (A2) 

For real s S (4m^,oo) one has 

c;,(s±iO) = — v.p. r ^:^iM±ij-^(,). (A3) 

The physical value of these functions corresponds then to the +iO prescription. In mathematical language, each 
function Qi is Hilbert transform of the corresponding function J-i . 

The discontinuities J-i{x) needed for the construction of Sec. |TV]are expressed in terms of the following two loga- 
rithmic functions 

L(.)=logi-44> (A4) 

M(s) = -2 log 1 + ^ +log (A5) 

\ s s / s 

where A(s) = A,,t(s) and cr(s) = (Tt^{s) were defined in ([T6l) and (|T8|) . respectively; and 

m± — rriri ± TOtt- (A6) 

For the logarithms we place the branch cut along the negative real axis and Imlogz G (— 7r,7r). The function L{s) is 

then real on the physical region. This is, however, not true for M(s), which is real only for s > (this corresponds 
to the TyTT — >■ TTTT scattering region). We have 

Tiix)^a{x), (A7) 

T2ix) = Lix), (A8) 

-^3(x) = ^, (A9) 
xa(x) 

•^4(a.)-cT(.)^§^, (AlO) 

= (All) 



and the numbers of subtraction taken for them in (IA1|) read ki — k2 ~ 1 and = ki = kc, = 0. Let us note that in 
these expressions the branch of the square root \^^^{s) is inessential. 

To find an analytic form of the Hilbert transform from its integral definition is a non-trivial task. For some functions 
we can use the roundabout way using the formula (jA2[) trying to find a function analytic in the complex plane except 
the branch cut on the interval (4mJ, oo) where it has the discontinuity equal to the value of the function J-i{s). All the 
functions satisfying this requirement differ just by polynomials that can be restricted by the UV and IR asymptotics 
of the integrals depending on the number of subtractions. 

It is easy to find that 

disc ( cr(s) log ^^4—^1 = - 4TO^)7rcr(s) = e{s - Aml)TrTi{s). (A12) 
Since we have defined Gi with one subtraction, we conclude that 

Cyi(s) = i(2 + a(.)log44TT) = IS^^-^^)- (^13) 

TT \ (T(s) + 1 / 

As we have seen, Gi{s) is connected with the 0{p'^) unitarity part of the amplitude corresponding to the single two- 
pion rescattering in the final state, it is therefore no surprise that this result restores the one-loop function J(s) which 
is a once subtracted scalar bubble with mass m^. 
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t^2(.) = ^log^^^, (AM) 
27r a[s) + 1 



Similarly, one arrives at 



Gsis)^- --log )', [ log^ ^ 1 + ^ ■ (A15) 

Note that these three functions have appeared also in the two-loop pion scattering computation (67| . 

For functions Giis) and ^5(5), this roundabout way does not work and we have to employ their integral repre- 
sentations (jAll) and (IA3|) . We have two possibilities how to compute them: either by integrating these expressions 
numerically or by means of the construction of analytic approximations which is described in Appendix B. One should 
note that having these two functions only in either of these approximate forms does not mean that they are worse 
than the others — since they depend only on two masses TOtt and m,, and on the variable s, once we fix the masses 
we can tabulate them. 

From these five functions we can obtain all the other functions appearing in the dispersive representation by using 
the following properties of the general Hilbert transform with n subtractions 



s 



TT 



da; /C(x) 



In order to simplify the following relations, we define 



n ..in 

s — X 



1. The formula for raising the number of subtraction is 



n+l n I (s)\ 

H (s)=H(s)-s"limJ^ . (A18) 



In terms of n-th derivation of (s) with respect to s, we can write this expression also as 

n+l n s" 9" " 

n is)^nis)-^g^nm. (ai9) 



2. Hilbert transform H(a) (*) of function 

= ^ (A20) 

(where A(s) is from (fT6|) equal to (s — m^)(s — rn^)) using the knowledge of the Hilbert transform H"'{s) of 
function /C(s) reads 

" , , h {s,m+,n)- h {s,m_,n) 



3. Obviously, we can lower the number of subtraction for /C(a) (s) 



"-1 ml h js, ml, n)-mlh {s, m?_,n) 

T-L(\){s) = 1-^ 2\ {All) 

^ ' s[m-^ — m_) 



and 
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4. In the case of same-mass particles, there appear functions 



Their Hilbert transform is 



n^)(') + y-^ (A25) 



5. Again by lowering the number of subtractions we arrive at 



nUs)= '^'^'f-''\ (A26) 
'V^^^)(,) = !Ll:i±!il!i_i). (A27) 



Appendix B: Functions Q4 and Gs, analytical approximations 

In this appendix we discuss in more detail the functions G4{s) and ^^5(5) from the previous appendix, which we 
do not know analytically. We first find their relation to some explicitly known analytic functions, which allows us 
to express them through more simple dispersive integrals, and then we construct an analytical approximation to the 
latter. 

We start by defining the following functions for complex z 



m\ — z 



ar,T,{z) = ^ , (Bl) 

A / m?_ — z 

>^i/2{z) = {mi - z)ann{z) = ^ rn?_ - z^ m\ - z , (B2) 

a{z)= lim a^^{z) = Jl^^^ , (B3) 

mn->-m„ V Z 

where the square roots are defined as vT — 1 with the cut along the negative real axis. Therefore, (T,,^- (z) has a cut 
along (to?_, m^) with the boundary values on the real axis 



ik^.(x)|i/2 ml<x<ml (B4) 
|cT^^(a;)|V2 x>m\ 



and the discontinuity equal to 

disccr^^(a;) = e{x - mi)e{ml - a;)|cr,,^(.x)|i/2. (B5) 

The function a{z) has the similar properties obtained by replacing m_ — ^ and m+ — ^ 2to^ in the previous formulae. 
Analogously, 

r |Ai/2(x)|i/2 x<ml 
Ai/2(x + iO) = <^ _i|Ai/2(x)|i/2 ml<x<ml , (B6) 

discAi/2(a;) = -0{x - mi)e{m\ - a;)|Ai/2(x)|^/^ (B7) 
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With the help of these elementary functions we can construclF^ 

F,iz)^aiz)log(^^^^y (B8) 

Al/2(z) \ Vc^r,7r(2) - 1/ \^ 1 - — (7^^ (z) ^ J 

Taking the principal branch of the logarithm with the cut along (—00, 0) and — vr < Imlogz < tt, their discontinuities 
are equal to 

- disc Fi{x) = e{x - \m\)u{x), (Bll) 

TT 

- disc F^{x) = e{x - 4ml)L{x), (B12) 

TT 

- disc I{x) = e{x) , (B13) 

where L(s) was defined in (IA4[) anc(3 for a; > 

f |A(a;)|i/2 x<m'^_ 
K/2{x) = < i|A(a;)|i/2 ml_ < x < m\ . (B16) 
[ |A(x)|i/2 a;>m^ 

The functions 7^4,5(2;) from (jA10p - (jAlip for x > Ami. can now be expressed in a convenient form as 

Fj{x) ^-I{x-iO)discFj{x) j=4,5. (B17) 

TT 



Using the formula 

disc/(x)g(x) = /(x - iO) disc + g(x + iO) disc /(x), (B18) 

we get then for j = 4, 5 

- disc [F,{x)I{x)] = e{x - Aml)F,{x) + + iO) . (B19) 

To find the functions with the desired discontinuities J^j_5(x), it remains therefore to evaluate the dispersion integrals 

By using (IB19[) and "dispersive integrations by parts" , there can be easily proved thatP^ 

g,{z)^-{F,{z)I{z)-Kj{z)). (B21) 

TT 

This representation of the functions Qj{z) for j = A^b allows to construct a systematic analytical approximation of 
them based on the approximations of much simpler integrals Kj{z). 



^^Note the relation of the functions F4,{z) and ^5(2;) to the functions Qi{s) and Q2{s) from l|A13|l and l)A14ll . respectively. 

Note that contrary to Aiy2(^)i the values of \\/2{^) are not x — iO boundary values of any analytic function in the cut complex plane. 
We can however write 

-^i/2(^) = ^1/2 (2: + sign(a; — ml)iO) = \i/2 + sign(x — m^)iO) , (B14) 
where with the above choice of the cut of the square root, the function 



-^l/2{-2) = iy-z ~ mi^m^ - 2 (B15) 

is analytic with a cut along {—oo,'m?_) U (m^,oo). 

It may seem that as a consequence of I IB19II this relation holds only up to a polynomial which does not contribute to the discontinuity. 
However, taking carefully into account the necessary number of subtractions for each term on both sides of the relations, there can be 
proved that such a polynomial is in fact absent. 
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The integrals Kj{z) have the following structure 



dx 1 „ , , , , /■^'"x drjc 



^ ~ ^ ^ fM^))+ ^ . g.-(Kx)l), (B22) 

where the functions /j(cr) are given by (|BTT|) and ((BT2|) such that 

(B23) 

and the functions (/j (a) for real a > are defined as 

5,(a) = /,(ia). (B24) 

The key idea is now to find appropriate series expansions of the functions fj{(j) and gj{<y) in the variable a and 
integrate then term by term. This can be done analytically and as we will see in what follows and the results can be 
partially summed. 

In order to perform this task, we have to distinguish three different regions, namely, I: a; > 4m^, where we have 
< (T < 1 and fj{cr) have a convergent Taylor expansion in the domain \a\ < 1; II: 2m^ < x < 4toJ, where the 
function a{x + iO) = i\a{x + iO)| with \a\ < 1 and .g^dcrl) have a convergent Taylor expansion in powers of \a\ and 
finally III: < x < 2m^, where again a{x + iO) = i\a{x + iO)| but now with \a\ > 1 and 5j(|cr|) can be expanded in a 
convergent Taylor expansion in the variable l/|cr|. 

Furthermore, we split the region I into three subregions Ia,Ib,Ic corresponding to the different explicit form of the 
function \i/2ix) for x < m?_, x £ {m^, '^+) f^^id x > respectively (see (|B16[) ). In what follows we therefore write 
Kj{z) as a sum of the contributions of five regions 

Kj{z) = Kf\z) + Ky{z) + K]^{z) + iff (z) + iff (z), (B25) 

which are ordered according to increasing x. 

Let us illustrate the the general recipe using iff ^(z) on the region III. Here 

54(k|) = -2\<j\ arctan (^-^^ , (B26) 

55(|a|) = -2arctan2 (^i^^ (B27) 

and after the substitution u = l/\a{x + iO)| we have in this region (we use the shortcuts a± ~ (t(toj_), A = m^m^) 

(^) - fr^) , f , . (B28) 

Using the expansions 



fc=0 



we can write 



^-+^- h J" (^ + ^W-^)^(t + ,-2)(t + .z2) 

Note that the series for gj (i) converge absolutely, an in addition the partial sums of them have the integrable 
majorants l/2ulog[(l +u)/{l — u)] and l/21og^[(l + u)/(l — u)], respectively. Therefore, the sum and the integral 
are interchangeable. 

Apparently we have to calculate one extra integral for each fc, however, in fact all the integrals can be obtained 
easily from one such integrals. Indeed, let us define for complex w 

mI^,\w) = /' , , (B31) 
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so that 



Kf\z) = 1^ 
^ ^ ' A(T+cr_ 



-2 oo 



fe=0 



Then we can write 



1 5* 



c«=0 



using the "generating integral" M\\\(w\o) that is for complex a, |a| < 1, equal to 

dw 1 



(B32) 



(B33) 



Now, as a consequence, 

1 a*^ 1 



fc! 5a'' 1 — aw 



^"^/ a=0 i=o 



i! 9a* 



(0) 



III 



a=0 



(B35) 



This formula can be easily understood. m'^1{w) is by definition an analytic function with a cut along (0, 1) and the 
discontinuity 



discM/i'i^(t) =Q{t)Q(\-t) 



^{t + <J-^^){t + a-J) 
It is therefore determined uniquely up to a polynomial. Because 



(B36) 



(B37) 



we can identify M^^^{w) up to some polynomial P}^i{w) with w'^M^jj (w), 

The unknown polynomial P}ii{u!) can be fixed by imposing the requirement of the appropriate asymptotics for 
w — )• DO, where Mj^^{w) — ^ 0. Choosing 



(B38) 



a=0 



(B39) 



we subtract from w'^Mj^^^{w) just the singular (and finite) part at w — )■ oo which ensures the right asymptotics. 
By explicit integration we obtain the result 



tanh 



1 f ^w-i + al 



tanh 



-1 f VHzLV^z^ 



In conclusion, the result for this region is given by the formula 



fe=0 



i=0 



! do} 



a=0 



(B40) 



(B41) 



W = — <7{z) 
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which allows for systematic analytic calculation of successive approximations. 

The remaining regions can be treated in the same way as was described above, with the only difference that the 
series expansion of the integrand contains now also the odd powers of a and therefore two "generating integrals" M^"'' 
and N'f \ A = II, la, lb, Ic are needed instead of one. 

The above calculation of the integrals K-^{z)^ A ~ III, II, la, lb, Ic gives the result in the form of the series expansion 
in the variables w^{z), where 

n,^^\z) = -a{z)-^ = ^-±—^ (B42) 

w'\z) = -a{zf = ^'^^^~\ (B43) 
w^'^iz) = w'^'iz) = w''{z) = a{zf = ^ , (B44) 

in the general form 

oo 

Kf{z) = P^K(z))5] (ai^-^^Aff K(z)) (B45) 
where the first order polynomials P^iw) are 



pni(^) = 1±2!L , (B46) 

P"H=-i±^, (B47) 
pi^(u;) = p^^\w) = P^^{w) = (B48) 

and the coefficients a^j'^ and 6jf^ correspond to the expansions of the functions fj{<j) and gj{\/(T), j = 4,5, in the 

integrands within the region considered. Let us remind that M^"^ (w) can be obtained from the "generating function" 
in the general form 

H = wH4^^\w) -w^Y^ ^^^Mf{a-^) = w^Mf{w) - Qf{w), (B49) 



i=0 



where q[^\w) is a polynomial of order k, and similarly for the functions N^^^w). In their convergence region the 
series J2T=o '^i'^A^'^ ^^'^ Sfclo ^m-'a'^'' summed up reproducing the even and the odd parts of the functions 

/j(cr) and gj{l/a) or gj{a) in the integrands within the region A, e.g. for K^^{z) and for \w\ < 1 we get 

oo ^ 

"fc^A^'' ^^'^~r= ^^•^tan . (B50) 



k=0 " 

Note, however, that these convergence regions do not generally coincide, so that this partial summations can not 
be made simultaneously for all n. Nevertheless, for real z ^ x + iO, such a summation reproduces exactly either the 
imaginary or the real part of the corresponding integral Kj^{z). For instance, in the region III we have 

^f'(- + '^)- , I iOT , , g.-(K^ + iO)l), (B51) 



which givefl 



111. , , ^ . V ^ — 2 |(t(x)| arctan Icrfa:)! 

Im K\'' {x + iO - ^B{x)B{2mi - x) ' ^ ^' ' ^ ^' B52 

\^i'^(x) 



*Note that \^l'^{x) = X^^^{x) iov < x < 2M^. 
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The above summation of the terms is possible for |o'(a;)| > 1 and gives the following contribution to the whole integral 
(note that the remaining part which cannot be summed explicitly is real) 



^ 4t]ii(-(^)-')'=<^( - -i^)-') = 2y=^arctany3^(^-l-^^ mS^( - a{.y 

(B53) 

Thus, since 



the imaginary part of this "generating integral" is equal to 

ImM/°)( - <j{x + iO)-2) = ACT+a_-4|^7r6i(a;)6'(2rr4 - .t) — . (B55) 

a{xY - 1 Ai/2(a;) 

Here cr(a;)^ < and the imaginary part is fully reproduced. The same can be done for the other regions, where 
we need, however, |(t(x)| < 1. One finds that there again the imaginary parts (or the real parts) are completely 
reproduced. 

In the physical region we have < ^{x) < 1; therefore the summation for the regions I and II can be performed. 
The convergence of the resummed approximations is illustrated in Fig. 21 
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